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This  technical  guide  provides  an  explanation  of  the  techniques  used  by 
the  Laser  Branch,  Laser  Microwave  Division,  US  Army  Environmental 
Hygiene  Agency,  to  evaluate  non- laser  optical  sources.  Hazard  criteria 
and  spectral  data  reduction  techniques  are  explained.  Radiometric 
measurements  are  not  Included.  The  Laser  Microwave  Division  Spectral 
Weighting  Program  (LMDSWP--a  Fortran  V computer  program)  Is  presented  In 
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1.  BACKGROUND. 

r 

a.  The  Laser  Branch,  Laser  Microwave  Division,  US  Army  Environmental 
Hygiene  Agency  (USAEHA),  evaluates  hazards  from  not  only  lasers,  but  also 
broad-band  optical  sources.  Examples  of  such  sources  Include  searchlights. 
Infrared  missile  guidance  systems,  optical  beacons  and  hospital  ultraviolet 
lamps.  This  guide  explains  the  techniques  used  to  reduce  spectroradlometric 
measurements  of  such  sources.  To  evaluate  a broad-band  optical  source  such 
as  an  arc,  a lamp  or  an  array  of  lamps  as  are  found  In  hospital,  military  and 
industrial  equipment.  It  1s  necessary  to  detenmlne  the  spectral  distribution 
of  the  optical  radiation.  The  spectral  distribution  of  Interest  Is  that  of 
the  accessible  emission,  which  may  differ  from  that  of  the  open  arc  or  lamp 
due  to  filtration  by  a plastic  or  glass  window  or  by  other  optical  elements 
in  the  system.  The  final  hazard  analysis  of  an  optical  source  requires  the 
weighted  sum  of  several  spectroradlometric  parameters  to  estimate  total 
retinal  irradiance  and  biologically>we1ghted  corneal  and  skin  irradlance. 

How  these  weightt;d  sums  are  evaluated  Is  explained  In  detail  in  "The 
Evaluation  of  Optical  Radiation  Hazards"  by  D.  H.  Sllney  and  B.  C.  Freasler, 
Applied  Optics.  Volume  12,  pages  1-24,  January  1973  (reference  lOd). 
Uuraviolet  radiation  exposure  limits  are  provided  in  AR  40-46. 

b.  During  1976-1977,  a new  computer  program  --  The  Laser  Microwave 
Division  Spectral  Weighting  Program  (LMDSWP)  --  was  developed  by  the  Data 
Processing  and  Technical  Information  Services  Branch  of  USAEHA  and  is 
explained  in  this  guide  in  Appendix  A.  Appendix  B provides  a program  listing 

I and  a table  of  values  for  the  spectral  weighting  functions.  Appendix  C 

I provides  operating  instructions  for  using  the  program. 

I 2.  REQUIRED  RADIOMETRIC  DATA.  The  spectral  irradiance  E^  should  be  complete 

I from  200  nm  to  at  least  1400  nm.  For  fluorescent  lamps  and  many  arc  lamps, 

; little  infrared  radiation  beyond  1200  nm  exists  and  can  be  neglected  If 

instrument  capability  Is  limited  to  1200  nm.  The  spectral  Irradiance  E.  at 
) the  nearest  point  of  access  (usually  at  the  glas%  cover)  Is  of  Interest  In 

I assessing  potential  ultraviolet  hazards  to  the  skin  aTid  lye,  and  potential 

1 hazards  to  the  skin  from  the  entire  spectrum.  The  spectral  radiance  L^  Is  of 

I Interest  for  assessing  potential  hazards  to  the  retina  and  should  be  complete 

I from  400  to  1400  nm  (the  retinal  hazard  region).  Again,  the  actual 

) measurements  beyond  1200  nm  can  generally  be  neglected.  The  values  of  Lx  may 

I be  closely  estimated  from  Ex  values  and  source  dimensions.  Radiometric 

f quantities  are  defined  In  Appendix  D. 
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3.  MEASUREMENT  TECHNIQUES.  The  spectrum  of  an  open  arc  proca»s_  (e.g. , wejding 
arc),  an  arc  lamp,  a gas  discharge  lamp  or  a fluorescent  lamp  consists  of 

line  structure  plus  a continuum.  Significant  errors  can  be  made  in 
representing  the  spectrum  and  weighting  the  spectrum  against  a biological 
action  spectrum  if  the  fraction  of  energy  in  each  line  is  not  properly  added 
to  the  continuum.  The  first  panel  of  Figure  1 shows  a hypothetical  spectral 
recording  from  a spectroradiometer.  If  spectral  points  were  arbitrarily 
recorded  every  5 nm,  most  of  the  line-peak  recordings  would  be  missed.  The 
width  of  the  triangular  line  at  half  of  the  peak  is  called  the  band-width  of 
the  monochromator/spectroradiometer.  The  recommended  method  for  representing 
the  spectrum  in  tabular  form  is  to  provide  the  measured  spectral  irradiance 
[uW/(cm2«nm)]  of  the  continuum  at  regular  intervals  (typically  every  5 nm) 
and  then  list  separately  the  irradiance  (MW/ cm2)  ip  each  line.  The  latter 
values  are  determined  by  subtracting  the  continuum  spectral  irradiance  at  the 
spectral  line  from  the  peak  reading  and  multiplying  that  value  by  the 
bandwidth  of  the  monochromater  (typically  2 to  5 nm).  These  line  irradiances 
are  listed  separately. 

4.  SPECTRAL  HISTOGRAM.  For  graphical  illustrations,  the  continuum  and  line 
structure  is  recombined  by  the  computer  program  into  a histogram.  The 
spectral  divisions  of  the  histogram  most  accurately  present  the  spectral 
resolution  of  the  data.  If  the  spectrum  is  represented  in  5-nm  intervals, 
the  irradiance  of  each  line  is  divided  by  5 nm  and  added  to  the  continuum 
spectral  irradiance  value  in  that  5-nm  interval  in  which  the  emission  line  is 
located.  As  an  example,  we  may  wish  to  represent  the  spectrum  of  a mercury 
arc  by  having  points  at  300  nm,  305  nm,  310  nm,  315  nm,  etc.  The  5-nm  band 
centered  at  305  nm  (i.e.,  302.5  to  307.5  nm)  contains  the  303-nm  emission 
line  of  mercury;  likewise,  the  band  centered  at  315  nm  contains  the  313-nm 
emission  line.  Since  the  band  centered  at  310  nm  contains  no  emission  line 
of  mercury,  it  truly  represents  only  the  continuum.  Panel  3 of  Figure  1 
illustrates  a histogram  plot. 

5.  PROCESSING  OF  DATA.  Many  of  the  calculations  which  are  useful  in  hazard 
analysis  require  weighting  the  spectrum  against  a biological  action  spectrum 
(e.g.,  erythema  or  photokeri titis  action  spectrum,  the  photopic  response  of 
the  eye,  and  the  retinal -injury  action  spectrum).  The  LMDSWP  computer 
program  (Appendix  B)  was  developed  to  simplify  this  data  reduction. 

Normally,  there  is  little  error  Introduced  by  using  the  digitized  spectral 
irradiance  values  used  in  the  histogram  plot  for  this  process.  However,  if 
the  lamp  spectrum  is  changing  rapidly  at  the  same  location  where  the 
weighting  spectrum  undergoes  a rapid  change,  significant  errors  can  be 
introduced.  It  is,  therefore,  preferable  to  weight  the  line  values 
separately  when  one  type  of  lamp  is  routinely  evaluated.  The  mercury  lines 
(254,  293,  303,  313,  365,  405,  435  nm,  etc.)  are  found  not  only  in  mercury 
lamps  but  also  in  fluorescent  lamps.  Routine  hand-computing  techniques  and 
the  machine-computing  routine  separately  weight  two  spectra  and  then  add 
them  afterward. 


Technical  Gutde  - Hazard  Analysts  of  Broad-Band  Optical  Sources 


275  on  300  nm  325  nm  275nm  300  nm  325  nm 

INSTRUMENT  READING  CALIBRATION 


275  nm  300  nm  325 1 
IRRADIANCE 


Figure  1.  Hypothetical  Spectral  Data  Reduction.  The  Instrument  reading  Is 
► divided  by  the  calibration  factor  to  yield  the  spectral  Irradlance. 
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6.  ABSOLUTE  VERSUS  RELATIVE  SPECTRAL  IRRADIANCE  MEASUREMENTS.  Absolute 
spectral  irradiance  or  spectral  radiance  measurements  are  not  always 
essential.  Photometric-to-radiometric  conversion  factors  (lumen/watt  ratios) 
can  be  obtained  from  a relative  spectrum  of  the  lamp  taken  at  any  accessible 
location.  Provided  that  a luminance  and  Illuminance  measurement  are  made  at 
a point  of  interest,  the  absolute  spectroradiometric  values  can  be 
calculated.  This  latter  approach  is  often  preferable  since  illuminance 
(lm/cm2)  and  luminance  (cd/cm2)  measurements  can  be  made  rapidly  at  many 
accessible  points  of  interest.  Provided  that  the  spectrum  does  not  change 
from  one  point  to  another,  the  ultraviolet-radiation  hazard  and  retinal- 
injury  hazard  can  be  calculated  at  all  such  points.  Because  of  the  limits  of 
sensitivity  introduced  by  photocathode  noise  in  typical  spectroradiometers, 
it  is  often  not  possible  to  obtain  spectral  irradiance  values  at  some 
distance  from  the  source.  This  problem  is  particularly  acute  in  the  near 
infrared  beyond  1000  nm.  By  measuring  the  relative  spectral  irradiance  as 
near  to  the  source  as  possible,  this  problem  can  be  greatly  reduced.  The 
LMDSWP  computer  program  calculates  a 1 umen-to-watt  conversion  factor. 

7.  SPECTRAL  RADIANCE  MEASUREMENTS.  The  spectral  radiance  Lx  [yW/lcm^-nm)] 
is  the  spectral  irradiance  E x at  a point  of  measurement  divided  by  the  solid 
angle  of  the  source.  It  is  typical  to  measure  EX  from  one  bright  spot  on  the 
source.  This  can  be  accomplished  by  using  an  irradiance  meter  with  a very 
small  fiel d-of-view  or  by  occluding  all  but  a small  portion  of  the  source 
with  a mask.  The  solid  angle  of  the  source  then  becomes  the  exposed  area  of 
the  source  Aj  divided  by  the  square  of  the  distance  r from  this  area  to  the 

i nstrument: 


L = E/iJj  and  = As/r^ 


(1) 


8.  BIOLOGICAL  WEIGHTING  FUNCTIONS.  Several  potential  hazards  must  be 
evaluated  separately  for  different  spectral  bands  and  biological  weighting 
functions  must  be  employed  for  this  purpose. 

a.  Ultraviolet  Radiation  (200-400  nm).  The  potential  hazards  associated 
with  ultraviolet  radiation  exposure  of  the  eye  and  skin  must  be  considered  in 
each  of  two  spectral  regions.  These  regions  are  the  actinic  (or  UV-B  and 
UV-C)  region,  where  photokeri ti tis,  conjunctivitis  (as  with  "welder's  flash") 
and  erythema  (as  with  "sunburn")  are  the  health  hazards  associated  with  this 
form  of  radiation  exposure;  and  the  near  ultraviolet  (UV-A)  spectral  region, 
where  the  effects  are  not  well  known,  but  cataractogenesis  has  been 
suggested.  Cataractogenesis  may  also  result  from  UV-B. 

(1)  Actinic  UV.  Standards  for  exposure  to  the  eye  and  skin  developed  at 
USAEHA  and  now  recommended  by  both  the  American  Conference  of  Governmental 
Industrial  Hygienists  (ACGIH)  and  the  National  Institute  of  Occupational 
Safety  and  Health  (NIOSH)  of  the  US  Department  of  Health, Education,  and 
Welfare  have  generally  become  accepted  in  the  USA  — particularly  where 
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ocular  exposure  is  of  concern  (references  10a  - 10c).  The  limit  for  exposure 
is  based  upon  an  "envelope"  action  spectrum  for  photokeri ti ti s and  erythema. 
The  spectrum  from  the  source  at  the  nearest  accessible  point  is  weighted  by 
this  curve  (Figure  2 and  Table  1)  for  wavelengths  less  than  318  nm. 

The  weighting  formula  is:  Eeff  = (2) 

and  the  permissible  8-hour  limit  for  exposure  is  10“7  w/cm2  for  Eeff  (a  total 
corresponding  to  an  exposure  dose  of  3 mJ/cm2). 

(2)  Near-Ultraviolet  Radiation.  Criteria  for  limiting  personnel 
exposure  to  UV-A  radiation  (320-400  nm)  are  presently  based  upon  limited 
biological  data.  The  solar  irradiance  incident  upon  the  skin  of  an 
individual  out-of-doors  is  normally  1-4  mW/cm2.  The  level  of  1 mW/cm2  is 
often  used  as  a safe  exposure  limit  (references  10a  - lOd).  Summing  the 
spectral  irradiance,  Ei  , from  320  to  400  nm,  one  obtains  the  total  irradiance 
in  the  UV-A. 

b.  Visible  Radiation  (400-770  nm)  and  Near-Infrared  Radiation  (770-1400 

nm) . 


(1)  Blue  Light  Hazard.  The  exact  boundaries  for  light  (or  "visible 
radiation")  are  often  argued;  at  present,  the  International  Commission  on 
Illumination  (CIE)  sets  380-400  nm  to  760-780  nm  as  "visible."  However,  of 
principal  interest  in  most  USAEHA  special  studies,  is  the  effect  of  all 
radiations  from  400  to  1400  nm  that  reach  the  retina.  Except  for  small 
children  and  aphakics  (those  with  the  crystalline  lens  removed  by  cataract 
surgery),  so  little  UV-A  radiation  reaches  the  retina  that  retinal  exposure 
in  that  spectral  region  is  considered  insignificant.  Until  recently,  retinal 
injury  from  high-intensi ty  light  sources  was  thought  to  be  thermal  injury  to 
retinal  tissue.  In  the  past  few  years,  it  has  become  increasingly  evident 
that  a photic  effect  which  has  as  its  basis  a photochemical  (e.g., 
phototoxic)  reaction  is  responsible  for  threshold  light-induced  retinal 
injury  for  exposure  durations  exceeding  10  seconds  (references  lOd  - lOg).  The 
blue-light  wavelengths  near  440  nm  appear  to  be  by  far  the  most  hazardous. 
Although  laser  safety  standards  reflect  a photochemical  injury  hypothesis  for 
light  exposures  greater  than  10  seconds,  they  were  initially  based  on  very 
little  data  available  at  the  time  of  their  development  --  1973  (references  lOd 
and  lOh).  For  the  purpose  of  evaluating  noncoherent,  broad-band  sources,  it 
is  more  reasonrtble  to  develop  a standard  for  lamp  exposures  directly  from  the 
threshold  retinal  injury  data.  A blue-light  hazard  function,  B,  , was 
developed  at  USAEHA  from  the  data  of  Ham  and  is  given  in  Table  Z.  It  has 
since  been  proposed  as  a possible  future  TLV®  by  ACGIH. 


® TLV  - Threshold  Limit  Values  for  Chemical  Substances  and  Physical  Agents  in 
the  Workroom  Environment  with  Intended  Changes  for  1977. 
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Figure  2 ■ Recommended  Ultraviolet  Radiation  Exposure  Standard* 

This  figure  was  adapted  from  a figure  developed  and  published 
by  the  American  Conference  of  Governmental  Industrial  Hygienists 
In  "Threshold  Limit  Values  for  Chemical  Substances  and 
Physical  Agents  In  the  Workroom  Environment  with  Intended 
Changes  for  1977". 
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TABLE  1.  RELATIVE  SPECTRAL  EFFECTIVENESS  BY  WAVELENGTH 


Relative 

f 

Spectral 

Wavelength 

Effectiveness 

i 

(nm) 

(mJ/cm^)* 



200 

100 

0.03 

210 

40 

0.075 

220 

25 

0.12 

230 

16 

0.19 

240 

10 

0.30 

250 

7.0 

0.43 

254 

6.0 

0.5 

260 

4.6 

0.65 

270 

3.0 

1.0 

280 

3.4 

0.88 

290 

4.7 

0.64 

300 

10 

0.30 

305 

50 

0.06 

310 

200 

0.015 

315 

1000 

0.003 

* 1 mj/cm2  = 10-3  j/cni2 
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TABLE  2.  SPECTRAL  WEIGHING  FUNCTIONS  FOR  ASSESSING  RETINAL  HAZARDS  FROM 
BROAD-BAND  OPTICAL  SOURCES 


Wavelength 

(nm) 

Blue-Light 

Hazard  Function 

Bi 

Burn  Hazard 

Function 

^ 

400 

0.10 

1.0 

405 

0.20 

2.0 

410 

0.40 

4.0 

415 

0.80 

8.0 

420 

0.90 

9.0 

425 

0.95 

9.5 

430 

0.98 

9.8 

435 

1.0 

10.0 

440 

1.0 

10.0 

445 

0.97 

9.7 

450 

0.94 

9.4 

455 

0.90 

9.0 

460 

0.80 

8.0 

465 

0.70 

7.0 

470 

0.62 

6.2 

475 

0.55 

5.5 

480 

0.45 

4.5 

485 

0.40 

4.0 

490 

0.22 

2.2 

495 

0.16 

1.6 

500-600 

10[(450-M/50] 

1.0 

600-700 

0.001 

1.0 

700-1060 

0.001 

10[(^-700)/515] 

1060-1400 

0.001 

0.2 
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(2)  Retinal  Exposure.  To  consider  a retinal  injury  safety  standard,  we 
must  first  define  the  relation  between  retinal  levels  and  lamp  brightness  (or 
radiance).  The  retinal  ir radiance  Ep  is  related  to  the  source  radiance  Lj 
that  is  being  directly  viewed  and  is  independent  of  viewing  distance.  It  is 
also  influenced  by  the  transmission  t of  the  ocular  media  in  front  of  the 
retina  and  upon  the  pupil  diameter  de,  which  for  a bright  visible  source  is 
normally  less  than  3 mm.  The  relation  is: 

Er  = 0.27  Ls-T-dg^  (3) 

This  equation  may  be  used  to  calculate  the  retinal  irradiance  at  just  one 
wavelength  or  in  a narrow  wavelength  band  (e.g.,  blue  light),  or  it  may  be 
used  to  calculate  the  total  retinal  irradiance  from  400  to  1400  nm.  In  the 
latter  case,  the  spectral  radiance  distribution  Lx  must  be  weighted  against 
the  spectral  transmittance  of  the  ocular  media  to  obtain  an  average  or 
effective  transmittance  of  the  ocular  media,  (reference  lOh).  The  formula 
is: 


^eff  = 

(3)  Photometric  Values.  The  source  brightness  can  also  be  considered  in 
photometric  terms. 

The  luminance  Ly  of  the  source  is  found  by: 

Lv  = 683  £V;^-L;^*4X  (5) 

Using  this  formula  and  the  CIE  luminous  efficiency  function  Vx  which  has  a 
maximal  value  of  1.0  at  550  nm  where  the  radiometric-to-photometric 
conversion  factor  is  683  lumens/watts,  from  formula  (5)  the  luminance  is 
calculated.  Following  the  same  approach  for  total  illuminance  Ey  at  the 
point  where  E;^  was  measured  is: 

Ey  = 683  Jv^-Ej^-AX  (6) 

These  formulae  also  permit  one  to  calculate  the  luminous  efficacy  of 
radiation  from  the  lamps  in  lumens/watt.  If  tte  spectrum  is  weighted 
against  the  scotopic  (rod)  response  functions  Vx* , we  then  obtain  the 
scotopic  efficiency,  where  the  maximal  value  would  be  1.0  at  500  nm  (the  peak 
response  of  rods).  Although  photometric  quantities  are  not  normally  used 
solely  for  comprehensive  hazard  evaluation,  many  relatively  inexpensive 
measuring  instruments  may  be  used  as  a crosscheck  on  radiometrical ly  measured 
values.  Furthermore,  cosine-corrected  photometric  instruments  may  be  used  to 
cosine-correct  spectrally  measured  UV  data. 

(4)  Exposure  Limits.  Simplified  standard  limits  developed  at  USAEHA  and 
also  proposed  as  a future  TLV  by  ACGIH  for  broad-band  light  sources  provide 
the  following  limits  expressed  in  terms  of  the  source  radiance. 
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(a)  To  protect  against  retinal  thermal  Injury,  the  spectral  radiance  of 
the  lamp  weighted  against  the  function  (Table  2)  should  not  exceed: 

L(HAZ)  = 

where  L Is  In  W/(cm2*sr)  and  t Is  the  viewing  duration  (or  pulse  duration  If 
the  source  Is  pulsed)  limited  to  1 ys  to  10  seconds,  and  a Is  the  angular 
subtense  of  the  source  In  radians.  If  the  source  Is  oblong  (e.g.,  a tubular 
flash  lamp),  the  angle  refers  to  the  longest  dimension.  For  Instance,  at  a 
viewing  distance  r = 100  cm  from  a xenon  flash  lamp  of  length  * 50  cm  the 


approximate  viewing  angle  Is: 

a = 2[arctan 

(^/2r)] 

or: 

a'?  ^/r 

for  smal 1 a 

= 50/100  = 0.5  radian 


This  relationship  may  also  be  expressed  In  retinal  terms;  I.e.,  Irradlance 
and  Image  diameter.  The  USAEHA  retinal  hazard  function  Is  graphed  In  Figure 

3. 


(b)  To  protect  against  retinal  Injury  from  blue-light  exposure,  the 
Integrated  spectral  radiance  of  the  lamp  weighted  against  the  blue-light 
hazard  function  (Table  2)  should  not  exceed: 

Lp(HAZ)  = 100  J/(cm2»sr)  for  t <1C^  seconds 

, . <9) 

L(HAZ)  = 10  mW/(cm2*sr)  for  t >10^  seconds 

For  a source  radiance  L which  exceeds  2 mW/(cm2*sr)  In  the  blue  region,  the 
permissible  exposure  duration  t (max)  In  seconds  Is  simply: 

t (max)  = 100  J/(cm2«sr)  / L (blue)  (10) 

These  latter  limits  are  greater  than  maximum  permissible  exposure  limits  for 
440-nm  laser  radiation  (AR  40-46  and  ANSI  Z-136.1). 

c.  Infrared  Exposure  (770  nm  to  1 mm).  The  total  accessible  average 
Irradlance  In  the  Infrared  from  most  sources  should  be  kept  below  10  mW/cm2. 
This  value  Is  to  protect  against  either  retinal  Injury  or  cataractogenesis 
(reference  lOd).  The  IR-A  radiance  should  be  less  than  0.6/o. 

9.  /WBIENT  LIGHT  LEVELS.  It  Is  often  considered  useful  to  compare  the 
spectral  radiance  of  a lamp  to  the  same  spectral  radiance  of  natural  light 
sources.  Figure  4 Is  a plot  of  spectral  radiance  of  the  solar  disc  (rj  * 6.9 
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Figure  3,  USAEHA  Permissible  Retinal  Irradiance  for  Momentary  Viewing  of 
Extended  Sources  as  a Function  of  Retinal  Image  Size. 
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Summer  tun  near  Zenim. 


Morylond  and  Cohiornto 
f=n5  mW.cm 


Maryland  5:45  pm  (=105  mW.cm'^ 
MoTy\ond  neor  luntef  8:05  pm  5* 

— — ^ f=;2  mW.tm 


^ 10 


Californio  neor  lunief  8:15  pm  2-5 

6=18  mW.cm* 


350  400  450  500  550  600  650  700  750  800  850  900  950 
WAVELENGTH  (nm  ) 


Figure  4.  Spectral  Radiance  of  the  Summer  Sun  for  Two  Localities.  Note  the 

change  In  short-wavelength  spectral  Irradlance  as  the  sun  approaches 
sunset. 
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I 

, X 10‘^*sr)  and  of  an  average  blue-sky  condition.  To  calculate  the  spectral 

radiance  of  snow  at  noonday,  divide  the  uppermost  curves  by  50,000.  For 
example,  the  approximate  noonday  spectral  radiances  at  440  nm  are  (from 
' Figure  4): 

Sun:  1.3  W/(cm2*sr*nm) 

Snow:  2.6  x 10~5  w/(cm*sr*nm) 

Sky:  7 x 10"®  W/(cm2*sr*nm) 

Until  a thorough  understanding  of  chronic  exposure  hazards  to  the  eye  has 
developed,  one  should  be  concerned  about  ocular  exposures  to  the  levels 
exceeding  those  of  the  latter  two  sources  for  long  periods  of  time  (reference 
lOi). 
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APPENDIX  A 
PROGRAM  DESCRIPTION 

PROGRAM;  Laser  Microwave  Division  Spectral  Weighting  Program  (LMDSWP) 
LANGUAGE:  FORTRAN  (EXEC  8 run  stream) 

PURPOSE:  Calculate  spectral  irradiance  data  from  uncorrected  instrument 
readings,  weight  by  biological  and  filter  transmittance  functions,  and 
produce  tables  and  graphs. 

INPUT:  Description  and  control  cards 
Weighting  functions  data 
General  function  data 
Source  spectral  data 

OUTPUT : Print-out  of  certain  input  values 
Print-out  of  calculated  data 
Cal  comp  plots  of  certain  calculated  data 

SEQUENCE:  Each  set  of  data  cards  must  be  ordered  in  ascending  wave  length 
(card  columns  1-4) 


INPUT  DATA;  For  a column-by-column  description  of  the  Input  data  cards  see 
coding  forms  and  card  layout  sheets.  Of  critical  importance  to  LMDSWP  is 
proper  ordering  of  data  values  according  to  associated  wavelength.  The 
FORTRAN  name  and  printer  heading  for  wave  length  is  LAMBDA.  Comments  follow 
concerning  Input  data. 

a.  Biodeck.  The  unitless  functions  of  wave  length  set  forth  In  Table  1 
are  seldom  changed  and  are  used  routinely. 

b.  Calibration  Data:  For  each  of  seven  instruments  there  are  two 
'spectral  calibration  functions',  one  associated  with  a visible  light 
calibration  source  and  one  with  an  ultraviolet  source.  A spectral  calibration 
value  at  a given  wavelength  is  computed  by  dividing  the  instrument  reading 

by  the  source  irradiance. 
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PROGRAM  DESCRIPTION 


1.  Read  two  'project  description'  cards  for  use  in  the  printer  output 
heading.  The  cards  should  contain  the  project  number,  date  of  initiation, 
and  a brief  description. 

2.  Read  the  'calculation  control'  card,  which  contains  the  number  of  filters 
considered  (FORTRAN  name  NUMFIL  may  take  on  values  0,1, or  2),  the  form  of 
the  filter  transmission  functions  in  terms  of  the  number  of  data  fields  to 

be  read  per  filter  (FORTRAN  NOCOFT  = 1 for  precalculated,  2 for  raw  data  pairs), 
the  form  of  the  calibration  data  to  be  read  (integer  CALDAT  -0  for  raw  data 
pairs,  1 for  prenlculated,  2 for  program  to  generate  '1.0'  for  all  wave- 
lengths), the  weighting  calculation  to  be  performed  (integer  GENWEI  ■ 00-15; 
see  Table  2),  whether  a general  function  is  tobe  read(integer  GENFUN  ■ 1)  or 
the  program-initialized  function  of  'one'  for  all  wavelengths  is  to  be  used 
(GENFUN  = 0). 

3.  Read  the  fourth  'distance  factor'  card  for  adjustment  of  source  data  by 
the  inverse  - square  law.  There  is  a factor  (FORTRAN  DFU)  for  UV  wave- 
lenths  and  another  (DFV)  for  visible  wavelengths. 

The  first  four  control  cards  must  be  present  in  each  run.  The  data  yet  to 
be  described  are  optional,  and  depend  on  'calculation  control'  parameters. 

4.  Read  (unless  CALDAT  =2)  and  process  calibration  data. 

a.  If  data  are  precalculated,  multiply  by  DFU  for  wavelengths  less  than 
or  equal  to  300  nm  or  by  DFV  for  wavelengths  greater  than  300. 

b.  If  data  are  in  raw  form,  read  dividends,  then  divisors;  check  for 
wavelength  matches;  perform  divisions;  and  incorporate  distance  factors  as  in 
(a). 

5.  Read  biological  data  deck  (BIODEK):  S-,  U-,  V-,  V*-,  T-,  T-A-,  C-A-,  and 
A-LAMBDA  are  separated  into  wavelength-paired  arrays. 

6.  Read,  also  into  paired  arrays,  remaining  biological  data,  one  set  at  a 
time:  B-LAMBDA,  X-BAR  LAMBDA,  Y-BAR  LAMBDA,  Z-BAR  LAMBDA,  P-445,  P-535,  and 
P-575. 

7.  Read  filter  transmission  data,  if  any,  and  process  if  necessary. 

a.  If  data  are  in  precomputed  form,  read  wavelength  and  transnlsil on  values. 

b.  If  data  are  in  raw  form,  read  wavelength,  dividend,  and  divisor.  Perform 
division  and  store  result  as  in  (a).  For  a two-filter  case,  filter  two  is 
considered  first. 
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8.  The  general  function  array  FOFX  was  Initialized  to  all  'ones'  at  the 
program  beginning.  If  calculation  control  GENFUN=1 , read  general  function 
data  Into  FOFX. 

•'  9.  Read  data  for  the  source  under  study:  name  of  event,  solid  angle, 

and  spectral  data.  Disregard  data  with  wavelength  less  than  the  largest 
previously  read-wavelength.  Negate  wavelengths  of  cards  punched  with 
'PEAK';  this  distinction  Is  necessary  for  the  following  calculation. 

10.  Determine  the  DELTA  value  associated  with  each  wavelength,  for  use 
In  various  spectral  weightings. 

11.  A wavelength  Interval  for  source  readings  (or  Intermediately- 
calculated  values)  may  not  exactly  match  an  interval  for  biological  or 
calibration  data.  Therefore,  In  weighting  calculations  subroutine  INTERP 
may  be  called  to  use  source  deck  wavelengths  for  Interpolation  of  'bio'  or 
'cal'  values. 

12.  Perform  calculations  to  obtain:  Eifofx(X);  E,^;L>;E>»  F;^(1 ) ,Ex‘ F;^(2) , 
or  Ex‘Fa(1,2);  E,.(A);  G-LAMBDA;  Integrated  sums  of  source  irradiance 
multiplied  spectrally  by  filter  and  bio  functions;  and  percentages  of  total 
Irradiance  which  are  UV,  visible,  or  near  -IR  radiation. 

13.  Create  Cal  comp  plot  file  containing  header  record,  wavelengths,  and  data. 

14.  Print  the  calculated  data. 

15.  Read  and  process  the  next  set.  If  any,  of  spectral  data. 

16.  After  all  spectral  sets  have  been  processed,  create  a Cal  comp  plot 
tape  for  histogram  presentation  of  output. 
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TABLE  1.  Biodeck  Data 


No  calculation 


TABLE  2.  Values  of  GENWEI  - General  Weighting  Calculations. 
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APPENDIX  B 

COMPUTER  LISTING 


1*  C 

2*  C 

3«  C 

<»•  C 

S*  C 

6*  C 

7*  C 

8*  C 

9»  C 

10*  C 

11*  C 

12*  C 

\ 13*  C 

14*  C 

\ 15*  C 

[ 16*  C 

} 17*  C 

1 18*  C 

\ 19*  C 

; 20*  c 

I 21*  C 

* 22*  C 

1 23*  C 

1 24*  C 

25*  C 

26*  C 

27*  C 

j 28*  C 

I 29*  C 

30*  C 

31*  C 

32*  C 

33*  C 

34*  C 

' 35*  C 

36*  C 

37*  C 

38*  C 

39*  C 

, ..  "O*  c 

; 41*  C 

42*  C 

43*  C 

44*  C 

45*  C 

46*  C 


PROGRaM  NANE  LASER  HICROUAVE  DIVISION  SPECTRAL  UETOHTINO  PROGRAM 

JUSTIFICATION  REQUESTED  BT  LASER  HICROUAVE  DIVISION  2 JAN  1976 
BV  OIV  CHIEF  ROBERT  T.  WANGENANN 

JOB  NUMBER  UNASSIGNED  AS  OF  30  JUNE  1976 

COMPUTER  UNIVAC  1108-EXEC  8 

PROGRAMMER  ROBERT  LEE  SCHMITT 

DATE  COMPLETED  JUNE  1976 

CHANGES  NONE 

SUBROUTINES  INITIL 
INTERP 
SUM 
SUHl 
GLAHCO 
FILSUM 
SEQUEN 
PRTCON 
HEADIN 
BDREAO 
HPLOT 

ABSTRACT  PROGRAM  REDUCES  SPECTRAL  DATA  T /K  EN  FROM  VARIOUS  BROADBAND 
OPTICAL  SOURCES  PRODUCING  TABLES  OF  PERTINENT  INFORMATION 
AND  PLOT  GRAPHS. 

INPUT  RECORD  LAYOUT 

SOURCE  DESCRIPTION  CARD 
ONE  OF  TUO 

1-78  DESCRIOTION 

TUO  OF  TUO 

1-72  DESCRIPTION 


CALCULATION  CONTROL  CARD 
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h •• 


47* 

C 

48* 

c 

43* 

c 

5 

NUMBER  OF  FILTER  S 

50* 

c 

10 

NUMBER  OF  DATA  FIELOS  FOR  FILTER  ONE 

51* 

c 

15 

NUMBER  OF  DATA  FIELDS  FOR  FILTER  THO 

52* 

c 

20-65  BY 

5 FILTER  CALCULATIONS 

53* 

c 

70 

NUMBER  OF  FIELDS  FOR  CALIBRATION  DATA 

54* 

c 

74t75 

GENERAL  VEIGHTING  FUNCTION 

5 5* 

c 

56* 

c 

DISTANCE 

FACTOR  CARO 

57* 

c 

58* 

c 

1-9 

DISTANCE  FACTOR  FOR  U.TRAVIOLET  LIGHT 

53* 

c 

10-18 

DISTANCE  FACTOR  FOR  VISIBLE  UGHT 

eo* 

c 

SI* 

c 

62* 

c 

GENERAL  DATA 

CARD 

63* 

c 

64* 

c 

1-4 

WAVE  LENGTH 

6 5* 

c 

9-17 

DATA WHERE  DATA  IS  CALIBRATION 

66* 

c 

FACTOR#  8-LANBDA  • X-BAR  LAMBDA  t V>BAR 

67* 

c 

LAMBDA  t Z-BAR  UNBOA  • P445  LAMBDA  t 

68* 

c 

P535  LAMBDA  • P57 5 LAHBD  t COMPUTED 

63* 

c 

FILTER  DATA  t DATA  READINGS  ^ 

70* 

c 

71* 

c 

FILTER  DATA 

IN  RAW  FORM 

72* 

c 

73* 

c 

1-4 

WAVE  LENGTH 

74* 

c 

9-17 

DIVIDEND 

75* 

c 

25-33 

DIVISOR 

76* 

c 

77* 

c 

BIO- DECK 

78* 

c 

79* 

c 

1-4 

WAVE  LENGTH 

PO* 

c 

5-13 

S-LAMBD  A 

81* 

c 

14-22 

U-LAMBO  A 

82* 

c 

23-31 

V-LAMBD  A 

83* 

c 

32-40 

VP-LAM3  DA 

84* 

c 

41-49 

T-LAM33  A 

85* 

c 

50-58 

TA-LAMB  DA 

86* 

c 

59-67 

CA-LAMS  DA 

87* 

c 

68-76 

A-LAM3D  A 

88* 

c 

89* 

c 

EVENT  NAfE  CARD 

90* 

c 

91* 

c 

1-18 

EVENT  NAME 

92* 

c 

93* 

c 

SOURCE  SOLID 

ANGLE  CARD 

94* 

c 

95* 

c 

1-9 

OMEGA 

96* 

c 

97* 

c 

98* 

c 

END  CARD 

99* 

c 

00* 

c 

1-3 

END 

,01* 

c 

.02* 

c 

.03* 

c 
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PRINTtR  LAYOUT 


lOY*  C 

105*  C 

106*  C 

107*  C 

loe*  C 

103*  C 

110*  C 

111*  C 

112*  C 

113*  C 

113*  C 

115*  C 

116*  C 

117*  C 

iia*  C 

119*  C 

120*  C 

121*  C 

122*  C 

123*  C 

124*  C 

125*  C 

126*  C 

127*  C 

128*  C 

123*  C 

130*  C 

131*  C 

132*  C 

133* 

1 34* 

1 35* 

136* 

137* 

138* 

139*  C 

140* 

141* 

142* 

143* 

144* 

145* 

146* 

147* 

148* 

149*  C 

ISO*  C 

151*  C 

1S2* 

153*  C 

154* 

155*  C 

156* 

157* 

158* 

159* 

160* 


2- 6  HAVE  LENGTH 

9-17  CALZ8RATT0K  FA  a OR 

20-28  INSTRUNENT  READINGS 

31-39  ADJUSTED  INSTRUMENT  READINGS 

42-50  SPECTRAL  IRRAOIANCE  OF  SOURCE 

53-61  SPECTRAL  RADIANCE  OF  SOURCE 

64-72  3-NM  SPECTRAL-RETINAL  IRRAOIANCE 

75-83  7-HN  SPECTRAL- RETINAL  IRRAOIANCE 

86-94  GENERAL  NEIOHTINC  FUNCTION  RESULTS 

97-105  SPECTRAL  FILTER  TRANSMISSIONS  FOR  FILTER  Oft 

108-116  SPECTRAL  FILTER  TRANSMISSIONS  FOR  FILTER  TWO 

119-127  SPECTRAL  FILTER  TR ANSMIS SI CNS  FOB  BOTH  FILTERS 

SINGLE  HAVE  CALCULATIONS 

3- 96  RESULT  DESCRIPTION 

100-108  RESULT 

FILE  RSGSHP-PLOT  LAY  OUT 

1-9  ELANBDA 


BIOOECK  TABLES 

REAL  SLAHBDI340»2I  8 ACGIH  UV  HAZARD  ENVELOPE  FUNCTION 

REAL  ULAMBOf  340tZt  a 1936  C IE  UV  SKIN  ERYTHEMA  ACTION  SPECTRUM 

REAL  ALAHaOI340t2i  a ANSI  Z 13 6 LASER  lEIGHTINC-UV  HAZARD  FUNCTICN 

REAL  TLAMBDI340t2l  a OCULAR  MEDIA  TRANSMISSION 

REAL  TALAMBI340i2l  a ABSORPTION  IN  THE  RETINA 

REAL  CALANBf 340t2l  a RECIPROCAL  OF  ANSI -NEAR-INFRARED  RETINAL 

BURN  CORRECTION  FACTOR 

REAL  VLAMB0l340t2»  a CIE-197D  PHOTOPIC  VISIBILITY  FUNCTION 

REAL  VPLAMBI340.2I  a CIE-19TD  SCOTOPIC  VISIBILITY  FUNCTION 

REAL  aLAMaDt340t2)  3 ACGIH  BLUE -LIGHT  HAZARD  FUNCTION 

REAL  XBLAHBt340t2l  a CIE-1931  BLUE  CHROMATICITY  COORDINATE 

REAL  YBLAHB(34atZt  a CITE-1931  GREEN  CHROMATICITY  COORDINATE 

REAL  Z6LANB(340»2I  a CIC-1931  RED  CHTOMATICTTY  COORDINATE 

REAL  P445LBI34Qt2»  a DARTNALL  NOMORGRAM  ABSORPTION  COEFFICIENT  BLX 

REAL  P535LB(340t2)  a DARTNALL  NOMOGRAM  ABSORPTION  COEFFICIENT  GREEN 

REAL  P575LBI340t2l  a DARTNALL  NOMOGRAM  ABSORPTION  COEFFICIENT  FC-  -EC 

END  BIOOECK  TABLES 

REAL  OFV 

REAL  DFU 

REAL  0P(2I 

REAL  OMEGA 

REAL  EINSTRt340«21 

REAL  F0FXf34Ot2l 

REAL  CF( 340*21 

8-3 


a DISTANT  FACTOR  TO  ADJUST  VISIBLE  RAOIATION 
CALIBRATION  FACTOR 

a DISTANT  FACTOR  TO  ADJUST  UV  RADIATION 
CALIBRATION  FACTOR 

a PUPIL  SIZE  USED  IN  RETINA  CALCULATIONS 
a SOURCE  SOLID  ANGLE  FOR  SPECTRAL  READINGS 
a SPECTRAL  READINGS  OF  SOURCE  UNDER  STUDY 
a FUNCTION  TO  MODIFY  INSTRUMENT  READINGS 
a CALIBRATION  FACTOR 


161* 

REAL 

EirOFX(34nt2l 

a 

ADJUSTED  INSTRUMENT  READINGS 

lt2» 

REAL 

EL AHB0(340.2) 

a 

SPECTRAL  IRRAOIANCE  OF  SOURCE  UNDER  STUDY 

ie3* 

REAL 

LL AHBDf 340t2I 

a 

SPECTRA,  RADIANCE  CF  SOURCE  UNDER  STUDY 

164* 

REAL 

ERETL8l340t3l 

a 

SPECTRAL -RETINAL  IRRADIANCE  FOR  3 AND  7-NM  PIPIL 

lfc5» 

REAL 

EOELLB 

a 

TOTAL  SPECTRAL  IRRAOIANCE  OF  SOURCE 

1S6* 

REAL 

LDELLB 

a 

TOTAL  RADIANCE  OF  SOURCE 

167« 

REAL 

ACCIH 

a 

EFFECTIVE  UV  RADIATION  ACCORDING  TO  THE 

16B* 

C 

AC6IN  STANDARD  ACTION  SPECTRUN 

169* 

REAL 

CTE 

a 

EFFECTIVE  UV  RADIATION  ACCORDING  TO  THE 

170* 

C 

1936  CIE  UV  ERYTICMA  ACTION  SPECTRUN 

171» 

REAL 

ANSI 

a 

EFFECTIVE  UV  RADIATION  ACCORDING  TO  THE 

172* 

c 

ANSI-ZI36  LASER  1C  16 HT TNG  UV  HAZARD  FUCTION 

173* 

REAL 

BLUHAI 

a 

BLUE  LIGHT  HAZARD  FUCTION  WEIGHTED  AGAINST 

174* 

c 

SPECTRAL  IRRAOIANCE 

175» 

REAL 

X6AR 

a 

1931  BLUE  CHROMA TICITY  COORDINATES  WEIGHTED 

176* 

c 

AGAINST  SPECTRAL  IRRADIANCE 

177* 

REAL 

YBAR 

a 

1931  GREEN  CHROMATICITY  COORDINATES  WEIGHTED 

178* 

c 

AGAINST  SPECTRAL  IRRAOIANCE 

179* 

REAL 

ZBAR 

a 

1931  RED  CHROMATICITY  COORDINATES  WEIGHTED 

ieo» 

c 

AGAINST  SPECTRAL  IRRAOIANCE 

iei» 

REAL 

P445 

a 

DARTNELL  NOMOGRAM  ABSORPTION  COEFFICIENT  FOR 

ie2» 

c 

BLUE  WEIGHTED  AGAINST  SPECTRAL  IRRADIANCE 

183* 

REAL 

P'35 

a 

DARTNELL  NOMOGRAM  ABSORPTION  COEFFICIENT  FO- 

ie4» 

c 

GREEN  WEIGHTED  AGAINST  SPECTRAL  IRRAOTANCE 

185* 

REAL 

P5  75 

a 

DARTNELL  NOMOGRAM  ABSORPTION  COEFFICIENT  FOR 

186* 

c 

RED  WEIGHTED  AGAINST  SPECTRAL  IRRADIANCE 

187» 

REAL 

VE 

a 

RADIANT  EFFACACY  OF  RADIATION  FROM  LAMBOA-MIN 

ie8» 

c 

TO  LAMBDA-MAX 

189» 

REAL 

VIE 

a 

FRACTION  CIE  SCOTDPTC  RADIATION  FROM  LAMBOA-MIN 

19a» 

c 

TO  LAMBDA-MAX 

191* 

REAL 

TRANS 

a 

EFFECTIVE  TRANSMISSION  OF  OCULAR  MEDIA  FROM 

192* 

c 

LAMBOA-MIN  TO  LAM30A-MAX 

193* 

REAL 

TRANTX 

a 

EFFECTIVE  TRANSMISSION  OF  OCULAR  MEDIA 

194* 

c 

MULTIPLIED  BY  SPECTRAL  ABSORPTION  OF  OCULA- 

195* 

c 

MEDIA 

196* 

REAL 

EECA 

a 

ANSI  LASER  MPE  WnGHTTNG  FACTOR  FOR  VISIBLE 

197* 

c 

AND  INFRARED-A 

198« 

real 

PCTUV 

a 

PERCENT  OF  TOTAL  IRRADIANCE  WHICH  IS  UV 

19  9* 

c 

RADIATION 

230* 

REAL 

PCTVI 

a 

PERCENT  OF  TOTAL  IRRADIANCE  WHICH  IS  VISIBLE 

201* 

c 

RADIATION 

202* 

REAL 

OCTNIR 

a 

PERCENT  OF  TOTAL  IRRADIANCE  WHICH  IS  NEAR 

203* 

c 

INFRARED  RADIATION 

204* 

REAL 

ILLUM 

a 

ILLUMINANCE  IN  LAMENS  PER  SQUARE  CENTIMETER 

205* 

REAL 

LUMTN 

a 

LUMINANCE  IN  CANDELAS  PER  SQUARE  CENTIMETE- 

206» 

REAL 

DATA 

a 

INPUT  BUFFER  FOR  TABLE  DATA 

207* 

REAL 

DATA! 

a 

INPUT  BUFFER  FOR  TABLE  DATA 

'’08» 

REAL 

WAVE 

a 

WAVELENGTH  ASSOCIATED  WITH  DATA  AND  OATAl 

209* 

REAL 

FTH340.2  1 

a 

FILTER  ONE  DATA  TABLE 

210* 

REAL 

FT2«34nt2) 

a 

FILTER  TWO  DATA  TABLE 

211* 

REAL 

EFT«  340t4) 

aSPECTRAL  FILTER  TRANSMISSION  FOR  FILTER  ONE  AND 

212* 

c 

FILTER  TWO  AND  FOR  BOTH  FILTERS 

213* 

REAL 

F ILTERtSI 

a 

EFFECTIVE  UV  IRRAOIANCE  ACCORDING  TO  THE 

214» 

c 

THREE  ACTION  SPECTRA  THROUGH  EITHER  FILTER 

215* 

c 

OR  BOTH  FILTERS 

216* 

REAL 

GLAH3D(340t2t 

a 

GENERAL  WEIGHTING  TABLE 

217» 

REAL 

DELTA( 3401 

a 

TABLE  OF  WAVE  LENGTH  INCREMENTS  FOR  EINSTR  DATA 
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218* 

REAL 

PREV 

8 

DURING  INPUT  OF  1C  AO INGS .CONTA INS 

THE 

219» 

C 

VALUE  OF  THE  LAST  CARO  READ 

220» 

INTEBER 

LAN80AI3*2I 

221* 

INTEGER 

FLTCNTI3t2l 

222* 

INTEGER 

CAL 

223* 

INTEGER 

GENE  UN 

a 

INDICATES  WHETHER  GENERAL  FUNCTION 

DATA 

22«l* 

c 

IS  TO  9E  READ 

225* 

INTEGER 

GENUEI 

a 

THE  GENERAL  WEIGHTING  FUNCTION  TO 

CALCUATE 

22S* 

INTEGER 

FILCALf  9> 

a 

INDICATES  WHICH  FILTER  CALCULATION 

S ARE 

TO 

22T* 

c 

BE  PERFORMED 

228« 

INTEGER 

CALH0RI14)  a 

DESCRIPTION  OF  CALI3IATI0N  DATA 

229* 

INTEGER 

0ESCRPI28I 

a 

DESCRIPTION  OF  SOURCE  UNDER  STUDY 

23a« 

INTEGER 

NONEIL 

a 

NUMBER  OF  FILTER 

231« 

INTEGER 

NOCOFl  a NUMBER 

OF  DATA  FIELDS  FOR  FILTER  ONE 

232* 

INTEGER 

N0C0E2  a NUMBER 

OF  DATA  FIELDS  FOR  FILTER  TWO 

233* 

INTEGER 

GENERL 

23<»« 

INTEGER 

CALOAT 

a 

INDICATES  WHETHER  CALIBRATION  DATA 

HILL 

BE  I (TUT 

235« 

c 

AS  RAW  DATA  OR  CALCULATED  DATA 

236* 

INTEGER 

EVENTIII 

a 

EVENT  NUMBER  OF  ai  N DATA  READINGS 

237* 

INTEGER 

EVENT2I12I 

23«» 

INTEGER 

PAGE 

a 

PAGE  COUNT  FOR  LISTING 

239« 

INTEGER 

DATEI21 

a 

DATE  OF  THE  RUN 

2m« 

INTEGER 

ENOHK 

a 

INDICATE  WHETHER  DATA  SET  HAS  ENDEC  CORRECTLY 

2«tl* 

INTEGER 

PRTLAMI 16.31 

a 

USED  IN  LIST  HEADING  FOR  GENERAL  WEIGHTING 

2»2» 

c 

FUNCTION  COLUMN 

243* 

I NTE  GE  R 

LSTUAV 

a 

CONTAINS  LAST  HAVE  LENGTH  DURING  CALCULATIONS 

244« 

INTEGER 

PEAK 

a 

INPUT  THAT  INDICA^'ES  WH^THEa  VALUE 

IS  A 

245* 

c 

PEAK  VALUE 

246* 

INTEGER 

3 ANPAS(3I 

a 

BAND  PASS  WITH  USED  WHEN  CONPUTINO 

WITH 

247* 

c 

a 

PEAK  VALUE 

248* 

INTEGER 

FUAVE 

249* 

DATA  I8ANPASII)  »I=1.3»/3t5 

flQ/ 

250* 

DATA  ( (PRTLAHf I* Jl t J=1 

.31 

1/ 

251* 

1»  (NOT  • 

.*USEO»  *.* 

• t 

252* 

2»  (S-L* 

.*AMBOAJ*. * 

• • 

253* 

3*  lO-L* 

t*AMB0A1*t  * 

• f 

2 5 4* 

4*  (A-L* 

»*  AMBDAJ  *.  * 

^ t 

255* 

5»  IT-L* 

,*AMBOAT*.* 

• # 

256* 

6»  (T-A-* 

• *LAM3DA*. * 1 

• t 

257* 

7*  (C-A-* 

•*LAHBDA*.*) 

• f 

258* 

8*  (V-L* 

.•AHBOA)*.* 

• • 

259* 

9*  (V*-L* 

.*AMBDAJ*.* 

• • 

260* 

A*  (8-L* 

.*AMBOAI*. * 

• t 

261* 

9*  IX- BAR* 

.*  LAMBO*»*A» 

• f 

262* 

C* I Y- BAR* 

.*  LAMBO*»*A> 

• • 

263* 

0*(2-BAR* 

.*  LAMBO*.*AI 

• f 

264* 

E*  IP-* 

.*4451  *.* 

265* 

E*  (P-* 

.*5351  *.* 

• t 

266* 

G*  IP-* 

.*5751  *.* 

•/ 

267* 

DATA  II 

LAM3DAII.JI.U: 

1.21.1=1.31/ 

268* 

1 *S- 

LA*  . *MBOA*. 

269* 

2 *U- 

LA*  . *MBOA*. 

270* 

3 * A- 

LA*  . *MBDA*/ 

271* 

DATA  1 IFLTCNTlIt J»  t J:i 

.21 

.1=1.31  / 

272* 

1 * FILTE*  • *R  ONE 

*. 

273* 

2 * FILTE*  t *R  TWO 

*. 

274* 

3 *BOTH  F*  t *ILTERS*/ 

B-S 


275* 
276* 
277* 
278» 
279* 
280* 
281* 
282* 
283* 
2 8‘*» 
285* 
286* 
287* 
288» 
289* 
299« 
291* 
292* 
293* 
29‘»» 
29S» 
296* 
297* 
298* 
299« 
303* 
331« 
302» 
303* 
30<4» 
305* 
33o« 
307* 
30S« 
333* 
313* 
311* 
312* 
313* 
3m» 
315* 
316* 
317« 
318* 
319« 
320* 
321* 
322* 
323* 
324« 
325* 
326* 
327« 
328* 
329* 
330» 
331* 


2 


C 

C 

c 

c 


c 

c 

c 


c 

c 

c 


c 

c 

c 


c 

c 

c 

c 

c 


087*  rvfwr?  / 
i*rr»d:*' t*Ncr 
CONTINUC 
DP(1  I : .09 
0PC2  » = .49 
HAVf  z 200 


•FIOURC*.*  . <8  SDtUr***  F SPCC*«*rML 

*T'.*  C«»*M  FOir»»  •••  • •/ 

a .lat  sautHEO 

a .TCN  SQUARCO 


INITIALIZE  GENERAL  FUNCTION  TO  UNITY.  IF  JMOTHCR  FUNCTION  IS  DESIRED 
TT  WILL  BE  INPUT  BY  CAROS  AND  AFFECT  ONLY  THE  SPECIFIED  RANGES 


DO  5 I = 1.240 
FOFXCI.U  = WAVE 

ForxiI.2l  = 1 

HAVE  = NAVE  * S 
5 CONTINUE 


READ  SOURCE  DESCRIPTION 


READ  10.  (OESCRPdl  »Irl.l4  I 
READ  10. lOESCRPI II .1=19.281 
10  FORMAT  113*6. A2I 


READ  CALCULATION  CONTROL  CARO 

REAO  la.NUMriL.NOCOFl.NOCOFZ.FILCAL.  CALOAT  iC ENVEI. 6ENFUN 
30  F0RMATI15IE) 

REAO  THE  CALIBRATION  OAT* 

REAO  40.0FU.DFV 
40  FORMAT « 2F9. 2) 

ITCAL=CALD*T+1 

GO  TC  C4S.1 30.126) .ITCAt 

45 — calibration  IS  RAU.  CF  IS  TO  9E  CALCULATfS. 

130-  calibration  is  PRE-CALCULATEO. 

126-  CALIBRATION  IS  ALL  ONE'S.  PROGR* H GENERATES  TABLE. 

45  READ  35.CALHDR 
35  FCRMATCl 3A6.A2) 


C THE  CALIBRATION  INPUT  IS  THE  RAU  DATA  AND  (T  NUST  BE  CONPUTED 

DO  6C  I : 1.341 

Rr*0l5.50.ERR  = 70)  CFf I.ll .CF  II .2  I 
50  F0RM*T(r4.a.T9.E9.2T 

60  CONTINUE 
70  READIO.GQI  ENOMK 
BO  FCRMAT(*3T 

IFCENDMK  .EG.  'END')  GOTO  100 
PRINT  90 

90  FORMATI 'O' .’CALIBRATION  OAT*  NOT  ENDED  CORWaLY’l 
STOP 

100  00  110  J : 1.341 

RE*OI5.S0.ERR=12OT  HAVE. DATA 
IF  I IFIXIU*VE*0.01I  .EQ.  IFIXI CF CJ .1 140.01 M GOTO  9S 
PRINT  95. WAVE 

95  FORMATI '0  WAVE  LENGTH  '.FS.O.'  OF  SPECTRAL  IRRADIANCE  DOES*. 


I 

.1 

1 
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i.- 


332* 

333* 

339* 

335* 

96 

336* 

337* 

99 

338» 

339* 

390* 

98 

391* 

392* 

393* 

97 

399* 

110 

395* 

120 

396* 

397» 

121 

398* 

399* 

125 

350* 

351* 

352* 

353* 

C 

359* 

C 

GE 

3 5 5* 

126 

356* 

357* 

358* 

359* 

360* 

361* 

127 

362* 

363* 

C 

369* 

C 

365* 

C 

386* 

130 

367* 

368* 

369* 

190 

370* 

371* 

372* 

373* 

135 

379* 

150 

375* 

160 

376* 

377* 

378* 

379* 

170 

380* 

C 

381* 

C 

382* 

c 

38  3* 
389* 
385* 
386* 

180 

387* 

388* 

1 ’NOT  MATCH  HAVE  LENGTH  OF  READING  FOR  C ALX8  RA  TION*  t 

2 ’FACTOR. •! 

STOP 

TFtOATA  .NC.  Of  GOTO  08 
PRINT  99*UAVE 

FORNATC’O  DIVISION  BY  ZERO  IN  CALIBRATCN  FACTOR  SECTION  AT’t 
1 *NAVE  LENGTH  ’tFS.OI 

STOP 

IFIUAVE  .OT.  3001  GOTO  97 
CF(J.2I  = CFfJ.2)  / DATA  • OFU 
GOTO  110 

CFfJtZI  = CFIJ.ZI  / DATA  • DFV 
CONT INUE 

IFII  .EG.  J1  GOTO  125 
PRINT  121 

FORMATI’O  CALIBRATION  FACTOR  ERROR  - NUMBER  RAN  DATA  NOT  MATCHED*! 
STOP 

READ (0.801  ENDMK 

IFIENDMK  .EG.  •END'I  GOTO  170 

PRINT  90 

STOP 

NERATE  CF  TA3LE  FROM  200- 1 AOO  ( UA  V EL  ENG  TNS  >- -A  LL  VALUES  = 1. 

CONTINUE 

VAL=195 

DO  127  1=1.391 
CF(  I.H=VAL*5 
VAL=CFII.l t 
CFI  1.21=1.0 
CONTINUE 
60  TO  170 

READ  COMPUTED  CALIBRATION  DATA 

READ  35.CALHDR 
DO  ISO  I = 1.391 

READ(5.190.ERR=160»  ( CF II . JT . J=1 .2  I 
F0RMAT(F9.D.T9.E9.2I 
IFICFII.II  .GT.  3001  CCTO  135 
CF(I.2I  = CF(I.2I  • DFU 
GOTO  150 

CFd.ZI  = CFd.ZI  • OFV 
CONTINUE 

READ (0.80  > ENDMK 

IFIENDMK  .EG.  'END*!  GOTO  170 

PRINT  90 

STOP 

CONTINUE 
READ  BIODECK 
DO  190  I = 1.391 

READ  (5.ie0.ERR=200  I HAVE.  SLAMBDI  I.  21  .ULAMBDI  I.  21  .VLAHBDI  1.2)  . 
1VPLAM3IZ.2I  .TLAMB0(I.2).T  ALAHBd.Z  >•  CALAMBII  .2  I • AL  AMBD  ( 1 . 2 1 
FORMAT(r9.0.eE9.21 
SLAMBDd.ll  = NAVE 
ULAMBDII.ll  = NAVE 

B-7 


389* 
390* 
391* 
392* 
393* 
i99* 
39  5» 
396* 
397* 
398* 
399* 
900* 
401* 
4l02« 
it03* 
■tOH* 
«0S* 
<<06* 
HOT* 
408* 
409» 
4134 
411» 
412» 
413» 
414* 
415* 
4 16* 
417» 
418* 
419» 
42a» 

421« 
422* 
423* 
424* 
425* 
426* 
427* 
428* 
429* 
430* 
431* 
4 32* 
433* 
4 34* 
435* 
436* 
437* 
438* 
439* 
440* 
441* 
442* 
443* 
444* 
445* 


VL«NB0(I«1>  = M8VC 
VPLAHBITtll  = HAVE 
TLAnBOdill  = NAVE 
TALAHBd*!)  = UAVE 
CALAn3d*ll  = NAVE 
ALANOOd*!  > = HAVE 
190  CONTINUE 
200  REAOIO.ZlOi  ENDHK 
210  F0RNATIA3I 

iriENONK  .EQ.  *ENO*l  OOTO  230 
PMNT  220tl 

220  r0RNATI*a*t«6I0DCCK  HAS  NO  END  C«D  AT  INAGE  ••131 
STOP 

230  CONTINUE 
C 

c 

CALL  BDR rADtBLANBOl 
CALL  BDRCADIXBLANBI 
CALL  aOREAOirSLAHBI 
CALL  90RCADd3LANai 
CALL  30RCA0tP44SLBI 
CALL  0CRCAniP535Lai 
CALL  3DRtAC«P57'-LBI 
C 

C READ  FILTER  TRANSNISSION  E UNC  TIO  NC  SI  * IF  ANY 

C filter  TMO  data  is  read  FIRST. 

C 

I : NUPFIL  ♦ 1 
50TCC450.290.390I.I 
C 

C determine  number  of  data  fields  FOR  FILTER  ONE 

c 

290  G0TCI3OO.36OI.N0COF1 
C 

C ONE  DATA  FIELD  FOR  FILTER  DNE 

C 

300  DO  320  r = 1.341 

REA0(5.310.ERRr330MFTld.Kl.K=l  .2  > 

310  FORMAT(F4.0.T9.E9.21 
320  CONTINUE 
330  READ (0.3401  ENDMK 
340  FCRMATCA3» 

IFCENOMK  .FQ.  'END*)  GOTO  450 
345  PRINT  350 

350  rORMAT(*0*. ’FILTER  DATA  HAS  NOT  ENDED  CORRCTLY*) 
STOP 
C 

C TUC  DATA  FIELDS  FOR  FILTER  ONE 

C 

360  00  380  I = 1.341 

READ  (5.370.ERR=330»  FTl I I.U .DATA .OAT A1 
370  FORMAT eF4.0. T9. 2E9. 21 
FTK  I.?>  = DATA  / DATAl 
380  CONTINUE 
GOTC  345 
C 

C DETERMINE  NUMBER  OF  DATA  FIELDS  FOR  FILTER  TWO 
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C 

<*47« 

390 

60T0f40Ot43altNOCOr3 

4«B* 

c 

BBS* 

c 

ONE  DATA  FIELD  FOR  FILTER  TUO 

450* 

c 

451« 

400 

DO  410  I = lf341 

453* 

REA0l5*310>ERR=430MFT3II«KI*K:n*3  ) 

4S]* 

410 

CONTINUE 

454« 

GOTO  345 

455* 

430 

READ(0t340l  ENOHK 

456* 

IFIENDMK  .EQ.  'END*!  GOTO  390 

457* 

GOTO  345 

4S8* 

c 

459* 

c 

TUO  DATA  FIELDS  FOR  FILTER  TUO 

480* 

c 

461* 

4 30 

DO  440  I = lt341 

463* 

READ  IS>  370  tERR=43a)  FT  311  >11 1 0 AT  A>  OK  TAl 

4S3* 

FT3II.31  = DATA  / DATAl 

464* 

440 

CONTINUE 

4E5* 

GOTO  345 

466* 

c 

4B7* 

c 

READ  GENERAL  FUNCTION.  INPUT  C0NTA3<S  THE  FIRST 

AND  LA  ST  UAVE  LENG  TH 

468* 

c 

AFFECTED  BY  FUNCTION 

469* 

c 

470* 

4 SO 

IF(6ENFUN  .EQ.  01  GOTO  530 

471» 

READ  15. 460.ERRr4e0l  FUAVE .LUAVE. DA  TA 

473* 

460 

F0RMATII4.T10.I4.T3O.E9.3 1 

473» 

c 

474* 

c 

OETE  RHINE  THE  TABLE  ELEMENT  NUMBERS  FOR  THE  UA  VE 

LENGTH  INTERVALS 

475» 

c 

476* 

J r (FUAVE  - 300»  / 5 ♦ 1 

477* 

K = (LUAVE  - 300)  / 5 ♦ 1 

478* 

c 

479* 

c 

INSERT  the  FUNCTION  INTO  THE  CORRECT  ELEMENTS  OF 

THE  GENERAL  FUNCTION 

480* 

c 

TABLE 

481* 

c 

483* 

DO  465  L = J.K 

483* 

F0FX(L.3)  = DATA 

484* 

465 

CONTINUE 

48S* 

GOTC  530 

486* 

480 

PRINT  510 

487* 

510 

F0RMAT(’0* .•INCORRECT  GENERAL  FUNCTION  CARD*  > 

488* 

STOP 

489* 

530 

CONTINUE 

4 9 3* 

CALL  lOENT 

491* 

c 

4934 

c 

READ  SPECTRAL  DATA  OF  SOURCE  UNDER  STUDY 

493* 

c 

494* 

c 

READ  THE  EVENT  NUMBER  OF  THF  DATA 

495* 

c 

498* 

530 

READ (5.540. EN0=9999)  EVENT 

497* 

5 40 

F0RMAT(3A6I 

498* 

c 

499* 

IF  (EVENTCII  .EG.  ’NEXT*!  GO  TO  3 

503* 

c 

READ  SOURCE  SOLID  ANGLE 

501* 

c 

503* 

READ  550. OMEGA 
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503 
' 0<» 
505 
505 

507 

508 

509 

510 

511 

512 

513 
51*« 
515 
51fi 

517 

518 

519 

520 

521 

522 

523 
529 
525 
5 26 
527 
529 

529 

530 

531 

532 

533 
539 

535 

536 

537 

538 

539 

590 

591 

592 
59  3 
599 

595 

596 

597 
5 9 3 
599 

550 

551 

552 

553 
559 

555 

556 

557 
559 
559 


C 

C 

C 


C 

C 

C 

C 


C 

c 

c 

c 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


95*0  THf  D*T*  REA0IN3S 
PRFV  : 0 

0 0 5 70  I = 1 #391 

555  95AOI5f  S60<rR9r5e0l  ( 5 INS  TR 1 1 1 Kt  •K  ^ I • PE  •( 

560  5 001**7159. 0tT3.C9.2. 720. A9) 

IF  PRESrN7  WAVE  LEN07H  IS  LESS  THAN  OR  ECU*.  TO  THE  PREVIOUS  WAVE 
LENOTH  disregard  IT  AND  READ  THE  NEXT  CARO 

iriPRFV  .LT.  EINSTRri.l»»  GOTO  561 
PRINT  562. «riNSTRII.Kl.K=1.2I 

562  FORMAT**  SOURCE  DATUM  DISREGARDED  - •.F9.0.E9.2I 
SOTC  555 

561  PREV  : EINSTRtl.H 

IF  THE  VALUE  IS  A PEAK  VALUE  THEN  NEGATE  T H"  WAVE  LENGTH  TO  INDICATE 
THAT  THE  VALUE  ASSOCIATED  WITH  THIS  WAVE  LENGTH  IS  A PEAK  VALUF 

IFCPEAK  .EG.  ‘PEAK*!  EINSTRII.11  = EINSTRII.ll  • t -1 » 

570  CONTINUE 
GOTO  600 

580  RE*0*0.59QI  ENDNK 

590  F0RMATIA3) 

IFCFNOMK  .FQ.  ’END**  GOTO  620 

600  PRINT  610 

Gin  format  * *0* . ’SPECTRAL  READINGS  WFRE  NOT  END  05  C0RRECTLY*1 
STOP 

620  CONTINUE 
NAX'TLMrl-l 

COMPUTE  THE  DELTA  ASSOCIATED  WITH  EACH  WAVE  LENGTH.  THIS  IS  DONE 
IN  THE  FOLLOWING  WAY! 

1)  DETERMINE  THE  OELFA  FOR  IHE  FIRST  NAVE  LENGTH.  IF  THE  WAVE 

LENGTH  IS  LESS  THAN  ZERO.  IT  IT  A PEAK  AND  DELTA  IS  SET  TO  CORRECT 
BANPA3  VALUE.  IF  IT  I S NOT  A PEAK  THEN  DELTA  EQUALS  ONE  HALF  OF 
THE  DIFFERENCE  BETWEEN  THE  FIRST  THO  CCNSECUTIVE  NON-PEAK  WAVE 
LENGTHS 

2»  COMPUTF  THE  REST  OF  THE  DELTAS.  PEAK  VALUES  GET  THE  CORRECT 
9ANPAS  VALUE.  FOR  NON-PEAK  VALUES.  MUST  KEEP  TRACK  OF  LAST 
NON-PEAK  WAVE  LENGTH  ANO  THE  WAVE  LENGTH  INTERVAL.  ALSO  MUST 
KNOW  IF  THE  LAST  VALUE  WAS  A PEAK  THAT  FELL  ON  A WAVE  INTERVAL 
WHICH  WILL  3E  TREATED  HAS  A NON-PEAK  FOR  THE  NEXT  NON-PEAK  DELTA 
VALUE.  SO  DELTA  FOR  NON  PEAK  VALUE  IS  T Ft  DIFFERENCE  BETWEEN 
THE  PRESENT  WAVE  LENGTH  ANO  THE  PREVIOUS  NON-PEAK  HAVE  LENGTH. 

IFJriMTRIl  .H.LT.O)  GOTO  623 
OIV  = 2 

WAVFl  : EINSTRII.ll 
LMN  = 3 

WAVr2  = EINSTR*2.1I 
UK  : 2 

619  DO  621  I : LMN.MAXELM 
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560* 

IFIHAVCl  .67.  0 .AND.  HAV  r2  .07.  0 1 6070  6S 

SSI* 

UAVn  r UAVC2 

SS2* 

HAVE2  = EINS7RfI»ll 

563* 

621 

C0N7INUE 

5S4* 

C 

SSS* 

c 

CONE  HERE  WHEN  7U0  C0NSECU7IVE  NON  PEAK  VALUES  ARE 

FOUND  AND  COMPUTE 

sse* 

c 

DEL7A  FOR  FIRS7  NON-PEAK  VALUE 

567* 

c 

588* 

622 

0EL7AI1I  =IUAVE2  - UAVEll  / OIV 

563* 

LS7UAV  = EINS7R(IJK-ltll  8 LAS7  NON-PEAK  NAVE  LEN67H 

570* 

IN7VL  = DEL7AI11  * DIV  8 INTERVAL  BETUFEN  2 CONSECUTIVE  NON-PEAKS 

571* 

6070  626 

572* 

c 

573* 

c 

FIRST  VALUE  IS  A PEAK  SO  DO  LOOP  UNTIL  A NON-PEAK 

VALUE  IS  FOUND 

574* 

c 

57  5* 

623 

DO  624  J = l.NAXELH 

576* 

DELTAIJI  = BANPASI3I 

577* 

IFIEINSTRf Jtll  .67.  -700)  DELTAIJI  = BANPASi2l 

578* 

IFIEINSTRI Jtll  .67.  -4001  DELTAIJI  = DANPASIll 

573* 

IFIEINSTRI J*l.ll  .67.  01  GOTO  625 

580* 

624 

CONTINUE 

581* 

6 25 

DIV  = 1 

582* 

UK  = J*2 

583* 

WAVE  1 = EINSTRI J*l.ll 

584* 

UAVr2  = EINSTRI J*2tll 

585* 

LNN  = J ♦ 3 

see* 

6070  619 

S87* 

c 

583* 

c 

PROGRAN  GOT  HERE  AFTER  IT  HAS  DETE  RHINE  DELTA  FOR 

THE  FIRST  NON-PEAK 

583* 

c 

VALUE  AND  ANT  PEAK  VALUE  THAT  WAS  BEFORE  IT 

530* 

c 

531* 

626 

DO  629  I = IJKtHAXELH 

532* 

IFIEINSTRir.il 1 627.631.578 

533* 

c 

594* 

c 

PRESENT  VALUE  IS  A PEAK.  IF  PEAK  I S ON  A NON-PEAK 

INTERVAL  then 

595* 

c 

PEAK  WAVE  BECOMES  LAST  WAVE  LENGTH 

596* 

c 

597* 

6 27 

IF  1 EINSTRII.il  - LSTWAV  .EO.  I NT  VL  1 LSTIRV  = EINSTRI!. 11 

598* 

DELTAIII  = BANPASfll 

599* 

IFIETNSTRII.il  .GT.  -7001  DELTAIII  = BANPASI21 

GOO* 

IFIEINSTRII.il  .GT.  -400)  DELTAIII  = 9ANPASIU 

601* 

80T0  629 

602* 

c 

G03* 

c 

PRESENT  VALUE  IS  NON-PEAK 

604* 

c 

605* 

6 28 

DELTAIII  = EINSTRII.l)  - LSTWAV 

GOG* 

LSTWAV  = EINSTRII.il 

607* 

INTVL  = DELTAIII 

608* 

629 

CONTINUE 

603* 

6 31 

CONTINUE 

610* 

c 

oil* 

c 

CALCULATION  SECTION  FOLLOWS 

612* 

c 

613* 

c 

614* 

c 

IT  IS  NOT  NECESSARY  THAT  THE  WAVE  LENGTH  INTERVAL 

OF  THE  SPECTRAL 

615* 

c 

READINGS  OR  ANY  OF  INTERMEDIATE  VALUES  CORRESPOND 

TO  THE  WAVE 

616* 

c 

LENGTH  INTERVAL  OF  THE  BIOLOGICAL  DATA  AND  THE  CALIBRATION  DATA. 
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617* 

C 

IT  IS  THEREEORE  NECESSARY  TH  EACH  CALCULATION  SECTION  TO 

NATCH  UP 

618* 

C 

BIO  AND  CAL  DATA  WITH  THE  CORRES POND INB  SPECTRAL  AND 

INTERMEDIATE 

619« 

c 

RESULTS  DATA.  IE  SPECTRAL  READING  FALL  BETWEEN  810 

DATA 

AND  CAL 

6rn» 

c 

DATA  THEN  THE  BIO  AND  CAL  DATA  WILL  NAVE  TO  BE  INTERPOLL ATEO  FOR 

621* 

c 

THE  HAVE  LENGTH  OF  SPECTRAL  DATAJ)0  TO  THE  FACT  THAT 

THE 

PROGRAM 

622* 

c 

NEGATES  THE  WAVE  LENGTH  OF  PEAK  READINGS  IT  IS  NECESSARY 

WHEN 

623* 

c 

COMPARING  WAVE  LENGTHS  TO  USE  THE  ABSOLUTE  VALUE  OF 

THE 

HAVE 

624» 

c 

LENGTH  FOR  SPECTRAL  READING  DATA 

625* 

c 

£26* 

c 

CALCULATE  ADJUSTED  INSTRUMENT  READINGS 

627* 

c 

628* 

K r 1 

629» 

DO  £40  I = ItMAXELM 

630* 

633 

IFfFOFXIKtl»  .EG.  ABS 1 FIN  STR(  I>1  >1  > GOTO  6S) 

631* 

IF(ABS(EINSTRIT.1I 1 .LT.  FOFXIKtlll  GOTO  638 

632* 

K r K ♦ 1 

S33* 

GOTO  £33 

634* 

£38 

EIFCFXIItll  = ABSIEINSTRf I«ll 1 

£35* 

CIF0FXII.2I  : EINSTRfT.21  • F0FX(K-lt2» 

636* 

GOTO  640 

637. 

£30 

EIFOFXd.ll  = ABSIEINSTRf  Ifll  1 

£38* 

EIFOFXIItZI  : EINSTRd.ZI  • F0FX«K.2» 

£39* 

K = K ♦ 1 

843» 

£40 

CONTINUE 

£41* 

c 

84Z» 

c 

CALCULATE  SPECTRAL  IRRADIANCE  OF  SOURCE  UNDER  STUDY 

£43* 

c 

£44* 

K = 1 

£45* 

DO  £90  I = l.MAXELM 

£4£* 

£45 

IFtEIFOFXL  ItU  .NE.  CFfKtlll  GOTO  £8  0 

£47* 

IF«CFIKt2»  .NE.  01  GOTO  S TO 

£48* 

£50 

PRINT  ££DtEIFOFX«Itll 

£49* 

££0 

FORMATC'D  ATTEMPT  TO  DIVIDE  A CALIBRATION  FACTOR  OF 

ZERO 

INTO  ADJU 

£50* 

ISTEP  INSTRUMENT  READING  AT  WAVE  LENGTH  •rf’S.OT 

£51* 

STOP 

£52* 

670 

ELAMBOfltll  = EIFOFXd.ll 

o53* 

ELAMBCd.21  = EIFOFXCI.ZI  / CFfK  .2  1 

£54* 

K = K ♦ 1 

£55* 

GOTO  £90 

£5£* 

£80 

IFCEIFOFXfl.ll  .LT.  CFfK. Ill  GOTO  £85 

£57* 

K r K ♦ 1 

££8* 

GOTO  645 

£59* 

£85 

CALL  INTERPI  EIFOFX.CFd.K  .CAL  FAC  I 

££3* 

IFCCALFAC  .EG.  01  GOTO  650 

££1* 

ELAM30d.ll  : EIFOFXd.ll 

£E2* 

ELAMB0d.2l  r EIFOFXd.21  / CALFAC 

££3* 

£90 

CONTINUE 

££4* 

c 

S£5* 

c 

CALCULATE  SPECTRAL  RADIANCE  OF  SOURCE  UNDER  STUDY 

£££• 

c 

££7* 

DO  700  I = l.MAXELM 

££8* 

LLAMBDd.il  = ELAMBDir.il 

££9* 

LLAMBOd.21  : ELAMBDd.21  / OMEGA 

873* 

700 

CONTINUE 

671* 

c 

£72* 

c 

CALCULATE  SPECTRAL-RETINAL  IRRADIANCE  FOR  3- MM  AND  7-HM 

PUPILS 

£73* 

c 
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67i»* 

8 = 1 

675* 

DO  750  I = ItNAXCLN 

676* 

705 

iriLLAHBDITtli  .NE.  TLAHBOfKtin  GOTO  720 

677* 

ERETLBIItlt  = LLANB0(I*1I 

678* 

J = 2 

679» 

DO  710  L = 1*2 

sea* 

ERETLBIItJ)  = .27  • LLAMBDri*2l  • ILAMBDIK*?!  • DPCLI 

681* 

J = J * 1 

582* 

710 

CONTINUE 

683* 

K = K * 1 

B8<** 

SOTO  750 

685* 

720 

IFILLANBOlItll  .LT.  TLANB0CK*1M  GOTO  730 

886* 

K = K * 1 

6a7« 

SOTO  705 

B88* 

7 30 

CALL  INTERP(LLAHBD*TLAHBO  *1*8 tTLRESt 

689* 

ERETLBlIfll  = LLANBDd.ll 

S90* 

J = 2 

691« 

DO  740  L = 1.2 

692* 

ERETL8(I.J»  = .27  • LLAHBD»I»2I  • TLRES 

* (T  TL  1 

693* 

J = J ♦ 1 

G9<»» 

740 

CONTINUE 

695* 

7 50 

CONTINUE 

698* 

C 

697* 

C 

CALCULATE  SPECTRAL  FILTER  TRANSMISSION 

FOR  ONE  FILTER 

698* 

C 

699* 

IFINUMFIL>  830*830.755 

700* 

755 

8 = 1 

701* 

DO  780  I = l.MAXELM 

702* 

760 

IFCELAMBOlI.ll  .NE.  FT1I8.11I  GOTO  765 

703* 

FTRES  = FT1I8.2I 

704* 

GOTO  775 

705* 

765 

IFIELAMBDII.il  .LT.  FTUM.HI  GOTO  7 70 

706* 

8 = 8*1 

707* 

SOTO  760 

708* 

770 

CALL  INTERPIELAMBD.FT1.T.8.FTRES  1 

709* 

775 

EFTII.ll  = ELAMBDII.il 

710* 

EFTII.2I  = ELAMBDII.2I  * FTRES 

711* 

8 = 8 *1 

712* 

7 80 

CONTINUE 

713* 

C 

714* 

c 

CALCULATE  SPECTRAL  FILTER  TRANSMISSION 

FOR  T UO  FILTER 

715* 

c 

716* 

GOTO (830.7e5I.NUHFIL 

717* 

785 

8 = 1 

718* 

J : 1 

719* 

DO  825  I = l.MAXELM 

720* 

7 90 

IFCELAMBDII.ll  .EG.  FTKJ.lll  GOTO  800 

721* 

IFCELAMBDII.ll  .LT.  FTlIJ.ll)  GOTO  795 

722* 

J = J *1 

723* 

GOTO  790 

724* 

795 

CALL  INTERPIELAMBD.FTl  ,1.  vl.FTlRESI 

725* 

GOTO  805 

726* 

8 00 

FTIRES  = FTHJ.2I 

727* 

805 

IFCELAMBDtI.ll  .EG.  FT2I8.1II  GOTO  815 

728* 

IFCELAMBDd.ll  .LT.  FT2IK  .111  GOTO  810 

729* 

K = 8 *1 

730* 

GOTO  805 
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7 31 

732 

733 

734 

735 

736 

737 

738 

739 

749 

741 

742 

743 
74  4 

745 

746 
74  7 
74  8 

74  9 

750 

751 

752 

75  3 
754 
75  5 

756 

757 
75? 

759 

760 

761 

762 

763 

764 

765 

766 

767 
769 

769 

770 

771 

772 

773 

774 

775 

776 

777 

778 

779 

780 

781 

782 

783 
7 84 
78  5 
786 
787- 


810  CALL  IN  TERP^CLA^•BDtFT2tIt  KtFT2RESl 
G070  820 

815  FT2RES  : rT2«Kf2l 

820  EF7tl.3I  : fLAMBDII.2»  • FT2RFS 

CFTIItRlr  rLAMB0(It2l  • FT2RES  • FTIRES 
K : K ♦ 1 
J : J ♦ 1 
825  CONTINUF 
830  CONTINUE 
C 

C DETERMINE  IF  A GENERAL  WEIGHTING  TABLE  IS  TO  BE  CALCULATED 

C 

3ENERL  = GFNWEI  ♦ 1 

GOTO ie6a«87n.e8nf89a.90a.9inf 920  t9  3nt940>95D  •9eO>970*9BO>99OtiaOOt 
11010 IfGENERL 

860  GOTO  1020  3 DO  NOT  CALCULATE  TABLE 

870  CALL  GLAMC0(GLAMB0.ELAM9DtSLAMB0*81O?O» 

880  CALL  GLAHC0(GLAMB0«ELAHBD •ULAMBDtSia2QI 
890  CALL  GLAMC0rGLAMBD.ELAM0DfALAMBO»*lO2O» 

900  CALL  GLAMC0IGLAMBDtELAMBDtTLAMB0.$1020» 

910  CALL  GLAMCOtGLAMB0.FLAM8D.TALAMB»*lD2OI 
920  CALL  GL AMCO t GLAMBOtELAMBD tCALAMB.«lQ20I 
930  CALL  GLAMCOJGLAMBD.ELAMBD tVLAMBD.*1020» 

940  CALL  GLAHCO (GLAHBDiELAHBO tVPLAMB tS 10201 
950  CALL  GLAMC0(CLAMBDtELAMBD«BLAMB0>S102OI 
960  CALL  GLAMCCIGLAMBDtELAM30tX3LAMBt$1020I 
970  CALL  GLAMCCCGLAMBD.ELAMBD tTBLAMB. *10201 
980  CALL  GLAMCOtGLAMBDtELAMOD t28LAMB.*1020l 
990  CALL  GLAMCO I GLAHEDtEL AM30 tP445LB t* 10201 
1000  CALL  GLAMGOIGLAMBDtELAMBD «P535L9»»1020> 

1010  CALL  GLAMCOIGLAMBD.ELAMBD tP575LB»*1020l 
1020  CONTINUE 

DC  1040  I r l.MAXELM 
IF  (EINSTRd.m  lOIOt  1040tl040 
1030  ELAMEDfl.l)  : EINSTRd.lJ 
LLAMPDd.il  : EINSTRd.il 

1040  CONTINUF 
C 

C !.».  JLATE  TOTAL  S»ECTRAL  IRRADIANCE  FOR  WAVE  LIMITS  USED  FOR 

C I Nd  hUMFNT  READINGS 

C 

ITFST  r 0 
GOTO  1045 

1041  DO  1141  ir=l.MAXELM 

1141  ELAMBDdI.2l:EFT(II.2l 
ITEST  : 1 

GOTO  1045 

1042  DO  1142  II=1.MAXELM 

1142  ELAMB0dI.2l:EFTdI.3» 

ITEST  r 2 

GOTO  1045 

1043  DO  1143  Iin.MAXELM 

1143  ELAMBDdI.21:EFTCII.4) 

ITEST  = 3 

1045  L r 1 

EDELL3  r SU  MU  E L AMBD  . OELT  A ,L  »M  AX  Q.  Ml 


c 

c 


CALCULATE  RADIANCE  OF  SOURCE  FRON  LANBOA-HTN  TO  LAHBOA-KAX 


I 


788* 
789* 
790* 
791* 
792» 
793* 
794* 
795* 
796* 
797* 
798* 
799* 
800* 
801* 
802* 
803* 
804* 
805* 
80&* 
807* 
808* 
809* 
810* 
811* 
812* 
813* 
314* 
315* 
’16* 
317* 
818* 
819* 
820* 
821* 
822* 
823* 
824* 
825* 
8ZB* 
327* 
828* 
8 29* 
830* 
831* 
332* 
833* 
834* 
83S* 
836* 
837* 
838* 
839* 
84  0* 
341* 
842* 
843* 
844* 


C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 

c 

c 

c 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


LOELLB  = EOELLB/OMEOA 

CALCULATE  EFFECTIVE  ULTRAVIOLET  RADIATION  C CORDING  TO  THE  ACOIH 
STANDARD  ACTION  SPECTRUM 

AC6IH  = SUM(CLAHB0*SLAMBD*0ELTA1 

CALCULATE  EFFECTIVE  ULTRAVIOLET  RADIATION  C CORDING  TO  THE  1936  CIE 
ULTRAVIOLET  ERYTHEMA  ACTION  SPECTRUM 

CIE  = SUHIELAHBDtULAMSDiDELTAt 

CALCULATE  EFFECTIVE  ULTRAVIOLET  RADIATION  C CORDING  TO  THE  ANST-Z13G 
LASER  WEIGHTING  UV  HAR2ARD  FUNCTION 

ANSI  = SUHIELAMBDt ALAMBO* DELTAI 

CALCULATE  BLUE  LIGHT  HAZARD  FUNCTION  WEI6HT0  AGAINST  SPECTRAL  IRRAOIANCE 
BLUHAZ  = SUMIELAMBDtBLAMBOtCELTA I 

BLUE fGREEN.RED  1931  CHROMATICITY  COORDINATES  WEIGHTED  AGAINST 
SPECTRAL  IRRAOIANCE 

XBAR  = SUMfELAHSOtXBLAMBt OELTAI 
YBAR  : SUMCELAMBD. YBLAMB. DELTAI 
ZBAR  = SUMIELAMBO.ZBLAMB. OELTAJ 

DARTNALL  NOMOGRAM  ABSORPTION  COEFFICIENT  FOR  B LUE*  GREE  N*  RED 
WEIGHTED  AGINST  SPECTRAL  IRRAOIANCE 

P445  : SUM«ELAMB0tP445L3t DELTA! 

P535  : SUMfELAMBDtP535LBt DELTAI 
P575  : SUMIELAMBD.P57SLB. DELTAI 

CALCULATE  RADIANT  EFFACACY  OF  RADIATION  FROM  LAMBDA-MIN  TO  LAMBDA-MAX 
VE  = 680  * SUM!ELAMBDtVLAMBD.DELTA  I / E0ELL3 

CALCULATE  FRACTION  CIE  SCCTOPIC  RADIATION  FROM  LAMBDA-MIN  TO  LAMBDA-MAX 
VIE  r SUMIELAMBDtVPLAMBtOELTA)  / EPELLB 

CALCULATE  EFFECTIVE  TRANSMISSION  OF  OCULAR  MEDIA  FROM  MIN  TO  MAX  LAMBDA 
TRANS  = SUMrELAMBD.TLAMBOtDELTA)  / EDELL8 

CALCULATE  EFFECTIVE  TRANSMISSION  OF  OCULAR  MFOTA  MULTIPLIED  BY 
SPECTRAL  ABSORPTION  OF  OCULAR  MEDIA 

TRANTX  = SUMIELAMBDfTALAMBtCELTA  I / EDELLB 

CALCULATE  ANSI  LASER  MPE  WEIGHTING  FACTOR  FT R VISIBLE  AND  INFRARED-A 
EECA  : SUMCELAMBD. CALAMBt DELTA!  / ECELLB 


I 
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S4S* 

34G* 

847* 

848* 

849* 

R5a« 

asi* 

852* 

853* 

854* 

855* 

856* 

857* 

ese* 

859* 
3Ba« 
861* 
862» 
863* 
a64» 
885* 
886* 
967* 
888* 
889* 
9 7 3* 
371* 
? 72* 
873* 
874* 
875* 
876* 
877* 
878* 
877* 
a 6Q* 
881* 
882* 
883* 
8 8.4* 
885* 
388* 
887* 
888* 
889* 
893* 
891* 
392* 
893* 
a94* 
895* 
898* 
897* 
898* 
999* 
903* 
901* 


ILLUH  = COriLB  • VC 
LUHIN  = ILLUH  / OHCO* 

C 

C CALCULATE  PERCENT  OF  TOTAL  IRRADIANCE  WHICH  IS  ULTRAVIOLET  RADIATION 

L = I 
1 = 0 
PCTUV  = 0 

IFIELAHBDIltlt  .GT.  4001  GOTO  1070 
DO  1050  I = IfHAXELH 
IFIELAHSDITtll  .GT.  4001  GOTO  lOGO 
1050  CONTINUE 

1060  IFJELAHBOII.il  .GT.  4001  1=1-1 

PCTUV  = 103  • SUHKELAHBD.OELTA.LtlT  / EDEU.B 
1070  CONTINUE 
C 

C CALCULATE  PERCENT  OF  TOTAL  IRRADIANCE  WHICH  IS  VISIBLE  RADIATION 

C 

IFJI  .EG.  HAXELHI  GOTO  1110 
HINELH  =1*1 

IFJHINELH  .LE.  II  HINELH  = 1 
'>CTVT  = 0 

IF(FLAH3D(T.1»  .GT.  7001  GOTO  1100 
DC  1030  I = HINELH. HAXELH 
IFJELAHBOJI.ll  .GT.  700)  GOTO  1090 
1080  CONTINUE 

1090  IFIELAHBDtI.il  .GT.  700)  1=1-1 

PCTVI  = 103  • SUHltELAHBO. DELTA. HINELH. 1)  / EOELLB 
HOC  CONTINUE 
C 

C CALCULATE  P'RCENT  OF  TOTAL  IRRADIANCE  WHICH  IS  NEAR  INFRARED 

C 

IFtI  .EG.  MAXELH)  GOTO  1110 
HINELH  =1*1 

IFtHINELH  .LC.  1)  HINELH  = 1 
PCTNIR  r 0 

PCTNIR=  103  . SUH1(ELAH30. DELTA. HINELH. HAXELH)  / EDELLB 
1110  IF  ( ITEST  .NE.  01  GOTO  1515 
C 

C CALCUATE  TOTAL  SPECTRAL  IRRADIANCE  WEIGHTED  SPECTRALLY  AGAINST 

C THE  TRANSHI5SI0N  OF  FILTER  ONE.FILTER  TWO  AND  BOTH  FILTERS 

C 

EEFl  = 0 
EEF2  = 0 
EEFIG2  = 0 

IFtNUHFILI  1330.1330.1170 
1170  GOTO (1180.1200) .NUHFIL 
1180  DO  1110  I = 1. HAXELH 

EEFl  r FEFI  ♦ EFTII.2)  • OELTA(I) 

1190  CONTINUE 
GOTO  1220 

1200  30  1210  I = 1. HAXELH 

EEF2  = EEF2  ♦ EFTtl.I)  * OELTAII) 

EEFIG2  = EFFIGI  ♦ EFT(I.4)  • DET-TA 
1210  CONTINUE 
3CTC  1180 
1220  CONTINUE 
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CALCULATC  THE  EFFECTZVC  ULTRAVTOUCT  IRRADZANCe  ACCORDZM  TO  ANY  OF 
THE  THREE  ACTZON  SFECTRAf  SLANIOt  UL  AMBDtALA  WOt  THROUW  EZTHER  FZLTER 
OR  BOTH  FZLTERS.  THERE  ARE  NZNE  POSSZBLE  OiANTZTZES  TO  CALCULATE 
ANO  THE  ONES  TO  BE  CALCULATED  ARE  SPECZFZED  ON  ZNPUT  BY  FZLCAL 


DO  1320  1 = 1*9 
CAL  = FZLCALIZI^l 
BOTO ll320tl230tl2R0tl2Sa>12Sa 
CALL  FZLSUBIFZLTERtCALtELAHBD 
CALL  FZLSUBIFZLTERtCALtELAMBO 
CALL  FZLSUBIFZLTERfCALtELAMBD 
CALL  FZLSUBIFZLTER>CAL*CLAHBO 
CALL  FrLSUBfFZLrERtCALtCLAMBO 
CALL  FZLSUBIFZLTER*CAL»ELAMBO 
CALL  FZLSUNIFZLTERiCAL>ELAMBD 
CALL  FZLSUNfFILTER»CAL*CLAHBO 
CALL  FZLSUNIFZLTERtCALiELAHBO 
CONTINUE 
CONTINUE 


1 12  TO  *1 2S0t  12  9 tl  300t  13 10 1*  CAL 
tFTlt  SLANB0»0ELTA>tl320l 
•FTltULANB0*0ELTAttl320> 

• FTl*  ILANBOtOfL  TA *<1320 1 
tFT2t  SLANBOtOEL  TA  »«1320l 
tFTZ*  ULANBOtOa  TA  *<13201 

• FT2*  ILANBOtOELTA  *<1320  I 

• FTl* FT2«SLANB0 (DELTA *<1320 1 
*FT1*  FT2*  ULAN  BD (DELTA  *<1320 ) 
(FTl*  FT2*  ALAN  BO (DELTA  *<1320  > 


ZOELTA  = ZnXIDELTAf  21 1 
URITE(1S*1331I  EVENT*NAXELN* ZOELTA 

1331  FORHATf 3A6*Z<*Z2I 

DO  1333  I = 1*HAXELH 

URirE(lS*1332l  ELANBOIZdl 

1332  FORNATIFS.Ot 

1333  CONTINUE 

DO  1335  I = 1*NAXELN 

URITE(1S*133<I  ELAH3DII*2I 
133<  FORHATieS.AI 

1335  CONTINUE 


PRINT  THE  TABLE  DATA NUST  DETERMINE  IF  TECRE  IS  FILTER  DATA  ANO 

IF  THERE  IS  A GENERAL  WEIGHTING  DATA  TO  BE  " RINTEO  ANO  IF  FILTER 
DATA  IS  TO  BE  PRINTED  WHETHER  IT  IS  ONE  OR  TWO  FILTERS 

CALL  DATEIHfDATEI 
PEAK  = 0 
PAGE  = 0 
LYNE:58 
LINMAX  =57 
K = 1 

IFIGENWEI  .GT.  01  K = 2 
K = K * NUHFIL  • 2 

G OTO II 340* 1370 *1400* lAlO* 1460*1490 )• K 

PRINT  TABLE  WITH  NO  GENERAL  WEIGHTING  ANO  NO  FILTER  DATA 


1340  K = 1 

00  1360  I = 1*NAXELH 
CALL  PRTCON 

CFRES  = SEQUENfEINSTR*I*CF*KI 

PRINT  13S0*EINSTRII*lt*CFRES*EINSTRI I* 21 *E IF OFXI 1*21 (ELAMBDI It  21  * 
1 LLAMBDII*2t*ERETLBII*2l  (ERETLBII*  31 


S59* 

1350 

FORMAino*  tT2>F5.a>Te>7llPC9.«>2X1  1 

9sa* 

1360 

CONTINUE 

961* 

80T0  ISIS 

9B2* 

C 

963* 

c 

PRINT  TABLE  UITH  GENERAL  WEIGHTING  DATAt  NO  FILTER  0i 

964* 

c 

965* 

1370 

K = 1 

9B6* 

DO  1390  I = ItNAXELH 

9C7» 

CALL  PRTCON 

9B8* 

CFRES  = SEGUENf EINSTRtI*CFtKI 

969* 

PRINT  1380.EINSTR(Itll  f CFRES  tEINSTRI  I.  2 1 tC  IF  OF  XI  It  21 

970* 

1 LLANB0IIt2ttlERETLBIItJI  tJ=2t3l  f0UHBDfIt2t 

971» 

1380 

FORMATI’O*  tF5.0tT8t8llPE9.4t2XM 

972» 

1390 

CONTINUE 

973* 

GOTO  1515 

974* 

C 

975* 

C 

PRINT  TABLE  UITH  ONE  FILTER  AND  NO  GENERAL  WEIGHTING 

97G» 

C 

977» 

1403 

K r 1 

978* 

00  1420  1 = ItPAXELM 

979* 

CALL  PRTCON 

980* 

CFRCS  : SEGUENIEINSTRtItCFtKI 

9ei» 

PRINT  IRlO.EINSTRdtlltCFREStEINSTRI  It  2 1 1 E IF  OF  XI  It  21 

982* 

1 LLAMBDI It  21 f lERE TLBI ItJI tJ-2t  31  tEFTIIt2I 

963* 

1413 

r ORPATfO*  f F5.0tT8t7(lPr9.4t2Xlt  T9  6t  1PE9.4  1 

984* 

1423 

CONTINUE 

965» 

GOTO  1515 

9 86* 

C 

38  7» 

C 

ORINT  TABLE  UITH  GENERAL  WEIGHTING  AND  ONE  FILTER 

988» 

C 

96  9* 

1430 

K = 1 

930* 

00  1453  I = ItMAXELM 

991* 

CALL  PRTCON 

992* 

CFRES  = SEQUENIEINSTRtTtCFtKI 

99  3» 

PRINT  144atEINSTR(It  1 1 tC'-REStEINSTRI  It  2 1 1 E 3F  OF  X 1 1 1 21 

994* 

1 LLAM30(It21t(ERErL3frtJI  ,jr2t  31  tGLAN8D(It21  tEFTIIt2 

39S» 

1440 

FORHATdO*  .r5.0tT8t9«lPE3.4t2XI1 

996* 

1450 

continue 

997» 

GCTO  1515 

993» 

C 

999» 

c 

ORINT  TABLE  UITH  TWO  FILTERS  AND  NO  GENERAL  UEIGHTTN 

inno* 

c 

3301* 

1460 

K r 1 

1302* 

OC  1430  T : ItPAXELM 

3 33  3* 

CALI  PRTCON 

13a4» 

CFRES  = SEGUENIEINSTRtItCFtKI 

13a5» 

PRINT  147atEINSTRIItl1tCPREStrTNSTRIIt21tEIF0FXIIt21 

1306* 

1 LLAM30IIt2l  f lERETLEI  It  J1  .Jr2t  31  tCEFTIItL1H.=2t41 

1307* 

1470 

FORMAT  1 ‘O’ tFE.0tT8t7llPE9.4t2XIt T9  6t  31 E9.4 t?  XI 1 

1038* 

1480 

CONTINUE 

1339* 

GOTO  1515 

1313* 

C 

13:i* 

C 

PRINT  TABLE  UITH  GENERAL  WEIGHTING  AND  BOTH  FILTERS 

1312* 

C 

1313* 

1490 

K = 1 

1314* 

DO  1510  I : ItMAXELM 

1315* 

CALL  PRTCON 
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1016* 

1017» 

1018* 

1019* 

lazo* 
lazi* 
lorz* 
1323» 
ia24* 
1325* 
1026* 
1027« 
1Q28» 
1329* 
1333* 
1331* 
1332* 
1333* 
1334* 
1335* 
1336* 
1337* 
1338* 
1339* 
134  3* 
1341* 
1342* 


CFRCS  = SCaUCNICINSTRtT*rr.KI 

PRTNT  ISOOtCINSTRIItlltCFRCStEINSTRI It  21 tC IFOFXI X* ?t  tE LAMBOt It  21 1 
1 LLAHaOIIt2ltlERETLBIItJ>  t.»=2t3}  fOUtHBOlItZI  tl  EFT!  It  LI  tL=2t4  I 
1530  F0RMATI*0* tF5.atT8tllllPE9.4t2XI  I 
1510  CONTINUE 
1515  PRINT  1520 
1520  F0HMATI*1*I 
PABE=PA0E*1 
PRINT  BOlOtOATEtPAGE 
PRINT  e020t0ESCRP 
PRINT  &025*EVENTt0ME6A 
6025  F0RHATIT2t'S0URCE  DESCRIPTION:  't 
A 3A6i*.  SOLID  ANGLE  = 'tES.!! 

IF  (ITEST  .EG.  01  GOTO  6130 
GOTO  l6131t6132t6133ltITEST 

6130  PRINT  6030 
GOTO  6040 

6131  PRINT  6031 
GOTO  6040 

6132  PRINT  6032 
GOTO  6040 

6133  PRINT  6033 

6010  F0RNATCT2t ‘LMD  SPECTRAL  WEIGHTING  PR  OGRAM*  .T  90  t • DA  TE  ’tZABt 
1 TllOt’PAGE  'tlTI 
C 

6023  FORHAT(*0«  t T2 1 1 3 A6 1 A2/ T7 1 13A6 1 A2  I 
C 


1043* 

6030  FORMATI //T45t*SUMMARIZATI0N  SHEET 

FO  R 

SOURCE 

1044* 

A / 

'tT2tl20l*-»n 

1045* 

6031  F0RMATI//T37t»SUMHARIZATI0N  SHEET 

FOR 

SOURCE 

1046* 

A WITH 

FILTER  0NE»t/tT2*12OC»-*Il 

1047* 

6032  F0RMATC//T37. •SUMMARIZATION  SHEET 

FOR 

SOURCE 

1348* 

A WITH 

FILTER  TWO*t/tT2tl20f*-*ll 

1049* 

6033  F0RMAT(//T3Qt*SUMMARIZATI0N  SHEET 

FOR 

SOURCE 

1350* 

A WITH 

FILTERS  ONE  AND  T WO  • t / 1 T2t  12  01  •- 

-*l  1 

1051* 

6040  PRINT 

leSOtEOELLS 

1052* 

PRINT 

1655tLDELLB 

1353* 

PRINT 

1660tACGIH 

1054* 

PRINT 

1670tCIE 

1355* 

PRINT 

1680tANSI 

1056* 

PRINT 

ISSOtBLUHAZ 

1057* 

PRINT 

ITOOtXBAR 

1058* 

PRINT 

ITlOtYBAR 

1359* 

PRINT 

1720tZBAR 

1360* 

PRINT 

1730tP445 

1361* 

PRINT 

1740tP53S 

1362* 

PRINT 

1750f P575 

1363* 

PRINT 

1760t VE 

13E4* 

PRINT 

ITTOtVIE 

1065* 

PRINT 

1780tTRANS 

1066* 

PRINT 

1790tTRANTX 

1067* 

PRINT 

IBOOtEECA 

1068* 

PRINT 

ISlOtPCTUV 

1069* 

PRINT 

1820tPCTVI 

1373* 

PRINT 

1830#PCTNIR 

1071* 

C 

1072* 

C PRINT 

FILTER  DATA  IF  PRESENT 
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1073* 
107«»» 
ia75» 
ia76* 
i077« 
1078* 
ia79» 
loaa* 
1081» 
1Q82* 
1083* 
iae4» 
1085* 
1086* 
1367* 
2 J8  8» 
I383* 
laao* 
1391* 
1392* 
2393* 
1394* 
1095* 
1396* 
1097* 
1098» 
1399» 
ii3a« 
1131« 
2102* 
2103* 
lia4» 
2205* 
210S* 
2137* 
2108* 
2239* 
1123* 
1111* 
1112* 
1113* 
1114* 
1115* 
1116* 
1117* 
1118* 
1119* 
1120* 
1121* 
1122* 
1123* 
1124* 
1125* 
1126* 
1127* 
1128* 
1129* 


IFfNUNFTL  .FG.  01  GOTO  1645 
PRINT  1840. FCFl 
G0T0fl540.153OI.NUMFrL 
1530  PRINT  1850. FEF2 

PRINT  1860.FCFIQ2 
1540  DO  1640  I r 1,3 

IFIFILCALdl  .CQ.  01  GOTO  1645 
N r FILCALCII 
GOTO  II 5 50.1 560,1570,1580, 

1550  PRINT  1870,L*MB0AI ll.FLTC 
GOTO  1640 

1560  PRINT  187n.LANB0*l2),FLT2 
GOTO  1640 

1570  PRINT  1870.LA>1BDAI3I.FLTC 
GOTO  1640 

1580  PRINT  1870,LAMBOACH.FLTC 
SOTO  1640 

1590  PRINT  187n,LAMB0AI2»,FLTC 
SOTO  1640 

1603  PRINT  1870,LAMBDAI3».FLTC 
GOTO  1640 

1610  PRINT  1870,LAMB0AI 1 ) .FLTC 
SOTO  1640 

1620  PRINT  1870.LAMBDAI2»,FLTC 
GOTO  1640 

1630  PRINT  1870.LAMBDAI3),FLTC 
1640  continue 

1645  PRINT  1880.ILI.UM 
PRINT  1890,LUMIN 
IF  INUMFIL  .EG.  0»  GOTO  1 
IF  (ITE5T  .EG.  01  GOTO  10 
GOTO  1 1042.10431 .ITEST 

1646  CALL  INITIL 
GOTO  530 

1650  FORMATdO  TOTAL  SPECTRAL 
1 FCR  THE  INSTRUMENT  REAOI 
1F,55  FORMATI’O  TOTAL  RADIANCE 
1 ,T130.1PE9.4I 

1660  FORMATI‘0  EFFECTIVE  ULTRA 
ISTANDARD  ACTION  SPECTRUM* 

1670  FORMATdO  EFFECTIVE  ULTRA 
ILTRAVIOLET  ERYTHEMA  ACTIO 
1680  F0RMATI*0  EFFECTIVE  ULTRA 
ILASER  WEIGHTING  UV  HA2AR0 
1690  FORMATI'D  BLUE  LIGHT  HAZA 
lACIANCE*  .T100,1PC9.4I 
1700  FORMATI'O  1931  BLUE  CHRPMATICTTY  C 
ITRAL  IRRAOIANCE*,T100,1PE9.41 
1710  FORMATI’O  1931  GREEN  CHROMATICITY 
ICTRAL  TRRA0rANCE*,Tin0,lPE9.4) 

1720  FORMATI’O  1931  RED  CHROMATICITY  CO 
IRAL  IRRADIANCE*,T100,1PE9,41 
1730  FORMATI’O  DARTNALL  NOMOGRAM  ASSORP 
ITED  AGAINST  SPECTRAL  IRRA OIANCE’ ,T 
1740  FORMAT! ’0  DARTNALL  NOMOGRAM  AQSORP 
ITED  A3AINST  SPECTRAL  IRRA 01 AN CE ’ .T 


1590,16  m .1  610,16  21  .16301  ,M 
NTfll  .FILTERtll 

NTI1I,F1LTERI21 

MTtll.F3LTERI3l 

MTI2I,FILTERI41 

MTI2I.FILTERI5I 

MTI2I,»'1LTERI6I 

NT!  31  ,FILTER«71 

MTI3l,Fn.TERIBl 

MTI31,FILTERI91 


646 

41 


TRRADIANCE  FOR  THE 
NGS’ ,T100 ,1PE9.4I 
OF  SOURCE  FROM  LA  16 

VIOLET  IR RADIANCE  A 
,T130,1PE9.4I 
VIOLET  IRRADIANCE  A 
N SPECTRUM’  ,T100,  IP 
VIOLET  IRRADIANCE  A 
FUNCTION’, T1 00, IPE 
•RD  FUNCTION  UEI0H1E 


WAVELENGTH  LIMITS  USED 

DA-MIN  TO  LAMBDA-MAX’ 

CC0R0IN6  TO  THE  ACGIH 

CCORDING  TO  1936  CIE  U 
E9.4I 

CCORDING  TO  ANSI-Z136 

9.41 

D AGAINST  SPECTRAL  IRR 


OOROINATTS  WEIGHTED  AGAINST  SPEC 
COORDINATES  WEIGHTED  AGAINST  SPE 
ORDINATES  WEIGHTED  AGAINST  SPECT 


TION  COEFF 
100,  1PE9.4 
TION  COEF*- 
100,1PE9.4 


ICIENT  FOR  BLUE  WEIGH 
I 

ICIENT  FOR  OREEN  WEIGH 
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113a*  1790  F0IWATI*0  OARTNALL  NOHOORIM  ABSORPTION  COCPFICICNT  FOR  RCO  WtTBH 

1131*  ITED  AGAINST  SPECTRAL  IRRA OIANCE* tT IQOt IPEB  I 

1132*  1760  F0RNATi«0  RADIANT  EFFICACY  ILUNENS/U ATT!  OF  RADIATION  FRON  LAMBDA- 

1133*  ININ  TO  LANB0A-NAX*iT100flPE9.*l 

113A*  1770  FORNATI'O  FACTION  CIE  SCOTDPIC  RADIATION  FROM  LANBDA-MIN  TO  LAMBDA 

1135*  l-NAX*tT10a*lPE9.A| 

1136*  1780  FORMATI’O  EFFECTIVE  TRANSMISSION  OF  OCULAR  MEDIA  FRON  LANBDA-MIN  T 

1137*  10  LAM80A-MAX**T1D0»1PE9.*I 

1138*  1790  FORMAT (*0  EFFECTIVE  TRANSMISSION  OF  OCULAR  MEOIA  MULTIPLIED  BY  SPE 

1139*  ICTRAL  ABSORPTION  OF  OCULAR  MEDIA ‘t  T1 00  tlPE  9.  «l 

IIRO*  1800  F0RNATr*O  ANSI  LASER  NPE  NEICHTINO  FACTOR  FOR  VISIBLE  AND  INFRARED 

1141*  1-A*f T10atlPE9.«l 

1142*  1810  F0RNATI*0  PERCENT  OF  TOTAL  IRRADIANCE  BETWEEN  LANBDA-MIN  AND  LANBD 

1143*  lA-NAX  WHICH  IS  ULTRAVIOLET  RAOIA  TION 't  T100*l  PC9.4I 

1144*  1820  FORNATC'O  PERCENT  OF  TOTAL  IRRADIANCE  BETWEEN  LAMBOA-MIN  AND  LAMBO 

1145*  lA-MAX  WHICH  IS  VISIBLE  RA DIATION 't TIOO tlPE 9l 41 

1146*  1830  F0RHATI«0  PERCENT  OF  TOTAL  IRRADIANCE  BETWEEN  LANBOA-NIN  AND  LAMBO 

1147*  lA-MAX  WHICH  IS  NEAR  INFRARED  RAOIATION*»Tiaa tlPCB.AI 

1148*  1840  F0RHAT(«0  TOTAL  SPECTRAL  IRRADIANCE  WEIGHTED  SPECTRALLY  AGAINST  TH 

1149*  IE  TRANSMISSION  OF  FILTER  ONE* tTlOO tlPC9.4 I 

USD*  1850  F0RMATt*0  TOTAL  SPECTRAL  IRRADIANCE  WEIOHTQ)  SPECTRALLY  AGAINST  TH 

1151*  IE  TRANSMISSION  OF  FILTER  TWO* t T1 DO *IPE9.4 I 

1152*  1860  F0RMATI»0  TOTAL  SPECTRAL  IRRADIANCE  WEIGHTED  SPECTRALLY  AGAINST  TH 

1153*  IE  TRANSMISSION  OF  BOTH  FI LTERS** T1 DO .1PE9. 41 

1154*  1870  FORHATMO  EFFECTIVE  UV  IRRADIANCE  ACCORDING  TO  •»2A4»  • ACTION  SPE 

1155*  ICTRA  THROUGH  *.  AGiTlDOtl PE9.4 I 

1156*  1880  FORHATI'O  ILLUMINANCE  IN  LUMENS  PER  SQUARE  C ENTIHETE R* tTlOD • 1PE9. 4 

1157*  II 

1158*  1890  F0RHATI*0  LUMINANCE  IN  CANDELAS  PER  SQUARE  C ENTIMETER* .TlOOt 1PE9.4 

1159*  II 

1160*  9999  END  FILE  15 

1161*  REWIND  15 

1162*  CALL  HPLOT 

1163*  PRINT  1900 


1164* 

1900 

FORMATI*!*.*  ALL  PROCESSING  COMPLETED*! 

1165* 

STOP 

1166* 

C 

1167* 

C 

• • 

1168* 

C 

# 

1169* 

C 

* END  OF  MAIN  PROGRAM 

1170* 

c 

1171* 

c 

• SUBROUTINES  FOLLOW 

1172* 

c 

0 

1173* 

c 

00 

1174* 

c 

1175* 

c 

00 

1176* 

c 

0 

1177* 

c 

* SUBROUTINE  IDENT 

1178* 

c 

• 

1179* 

c 

00 

1180* 

c 

1181* 

c 

1182* 

c 

PRINT  OUT  THE  INPUT  PARAMETERS  CONVERTING  70 

1183* 

c 

NARRATIVE  WHEN  POSSIBLE 

1184* 

c 

1185* 

c 

1186* 

c 

B-21 


1187* 

1188* 

1189* 

1X90* 

1191* 

1I92* 


SUBR OUT  INC 
DINENSION  r 
INTEGER  riL 
INTEGER  CAL 
INTEGER  GET 
DATA  I IFILT 


IDENT 

ILTA3lia*7l .CA 

TAB 

TAB 

AB 

ABIIT»JTI»JT=1 


LTABI  ?•  31  tSFTABdGt  Al 


JTI»JT=1.71 ,ITrl. 10  I 


1193* 

A ’NO  FIL 

►’TER  CA’» 

•LCUIAT 

•••ION  •• 

• • 

1194» 

I’E-L AMB 

t’DA*FIL’» 

•TER1«S 

•••LAMBDA’* 

••DELTAS 

1195* 

2’E-LAMB 

•’OA»FIL’» 

•TER1»U 

•••-L  AM  BO*. 

•A«DELr 

• •A 

1196* 

3’E-L  AMB 

.’DA«FIL’t 

’TCR1»A 

•••-LAMBO^. 

•A»DELT® 

• •A 

1197* 

4’E-L  AMB 

.•OA»FIL’. 

•TER2»S 

•••-LAMBD^. 

•A»OCLT* 

••A 

1198» 

S’E-L AMB 

t’DA*FlL’t 

’TER7»U 

•••-LAMBO^. 

•A*OELT® 

••A 

1199* 

6’E-L AMB 

.•DA»FIL’t 

•TER2»F 

•••-LAMBO^. 

•A^OELT 

• •A 

120D» 

7’E-L AMB 

t’DA»FIL’» 

•TER1»F 

•.•ILTER2^. 

••U-LAB* 

••BOA*DE 

•LTA 

1201« 

8’E-L  AMB 

t’DA«FIL’t 

•TER1*F 

•.’ILTER2 •• 

••U-LAB^ 

.•BDA*OE 

•LTA 

12Q2* 

9’ E-LAMB 

t’DA»FIL’. 

•TER1»F 

•.•ILTER2 •• 

••F-LAB* 

••BOA*OE 

•LTA 

1203* 

DATA  ( ICALTABIICtJCI* JC=1 

.3l*I(r=l*  2t 

/ 

1204* 

I’RAH  OA 

’.’TA 

9 

• • 

1205* 

2’CALCUL 

’t’ATED  O’. 

’ ATA 

•/ 

1206» 

DATA  1 ( 

GFTABIIGt JG 

I tJCrl. 

4)  .IG=1  *1  61  / 

1237» 

I’HEIGHT 

’.’ING  FN*. 

’ IN  PR 

•••06RAN  *9 

1208» 

2’ELAMBD 

’ t’ A»S-LA’t 

’MBOA 

* t * 

1209« 

3’ELAM3D 

’.•A»U-LA’. 

’MBDA 

* • * *9 

121'3« 

4’ELAM8D 

’ t’ A»A-LA’» 

’MBOA 

* 9 * * 9 

:2ii» 

S’ELAMBO 

’.•A»T-LA’. 

’MBOA 

* 9 * * 9 

1212* 

6’ELAMBD 

’ .’  A»TA-L’» 

’ AM30A 

9 g 9 9 y 

1213* 

7’ELAM3D 

’.’A»CA-L’» 

’ AM3DA 

* 9 * *9 

1214» 

9’ELAM3D 

• .’A’U-LA’. 

’MBOA 

* 9 * *9 

1215» 

9’ELAM3D 

’.•A»VP-L’t 

’ AM3DA 

* 9 * *9 

1216* 

A’ELAMBD 

’ t’A»D-LA’. 

’MBOA 

* 9 * * 9 

1217* 

O’ELAMCD 

• .’A^XBAR’t 

’-LAMBO 

*9**  *9 

:2i9» 

C’ EL  A MOD 

’ .’A  »Y3AR’f 

’-LAM3D 

♦9*A  *9 

1219’ 

D’ELAMSC 

’ .’A •ZBAR’t 

’-LAMB3 

’9*A  *9 

1220* 

E’ELAMBD 

• ,’A«P445’t 

•-LAM33 

*9*A  *9 

1221* 

F’ELAMBD 

’ t’ A»P535’t 

’-LAM90 

’9*A  *9 

1222* 

O’ELAMBO 

•.’A«P575’. 

’-L AM33 

♦9*A  •/ 

12r3»  6000  '•ORNATCIHI, 

6010  FORMATCTZO. 
1225*  ADATA  USED  F 

1226*  6011  F0RNAT«T20t 

1227»  6012  F0RHATIT23, 

1228*  A DATA  USEO 

1220*  6013  F0RMAT(T2O. 

1230*  6001  F0RMAT(T2O. 

1231*  6020  F0RMAT(T20. 

1232*  6021  F0RMATCT20. 

1233*  A»/» 

1234*  6022  F0RMAHT20. 

1235*  A 

1236»  6030  FORMAT! T20t 

1237»  A 

1238*  6031  F0RMAT(T20f 

1239*  AON  OATA»/l 

1240«  6040  F0RMAHT20. 

1241*  6050  rCRMAT(T20. 

1242*  6060  FORMAT  «T20 

1243*  A T20 


T5.*  DESCRIPTI 
’NCCCFl  : *.12 
OR  FILTER  ON"’ 
•NOCCFl  = ’.I? 
•N0C0F2  r ' ,Z? 
FOR  filter  TWO 
•N0C0F2  = '.r? 
’NUMBER  OF  FIL 
•CALDAT  : ’.I? 
•CALDAT  : ’.12 


ON  : ’t  13  A6 
t’  CALCULA 


t ’ RAW 
TTRS  USI 
f*  CAL] 
CAL] 


] fA7/.T  T 
ITED 

rA  USEO  *■ 
kTED 

rA  USED  *• 
: •tI2t/l 
ITION  I^P 
ITICN  IN® 


tl  3A6t  A?  < 


OR  FILIER  ONE’/I 


OR  FILTER  TUO’/» 


UT  IS  RA 
UT  IS  CA 


•CALOAT  = ••I2«’  NO  CA  LIBR  ATI  0 N INPU'.  t T 
•ONES’/ I 

•GENFUN  : •tT2t*  GENERAL  FUNCTION  IS  IN 
• PROGRAM’/T 

’GENFUN  = ’tI7t’  READ  INTO  COMPUTER  THE 

7A6» 

4A6I 

t ’ULTRAVTCLET  DISTANCE  : •tF9.2«/. 
t’VISIBLF  LIGHT  DISTANCE  = ’tFB^A 


U DATA’/) 
LCULATED  DATA 

ABLE  SET  TO  ’. 

COMPUTER’ 

GENERAL  FUNCTI 


1245* 
124G* 
1247* 
12<l8* 
1249* 
1250* 
1251* 
1252* 
1253* 
1254* 
1255* 
1256* 
1257* 
1258* 
1253* 
1260* 
1261* 
1262* 
1263* 
1264* 
1265* 
126  6* 
1267* 
1268* 
1269* 
1270* 
1271* 
1272* 
1273* 
1274* 
1275* 
1276* 
1277* 
1278* 
1279* 
1289* 
1281* 
1282* 
1283* 
1284* 
1285* 
1286* 
1287* 
1288* 
1289* 
1290* 
1291* 
1292* 
1293* 
1294* 
1295* 
1296* 
1297* 
1298* 
1299* 
1300* 


C 

C 

c 

C PRINT  DESCRIPTION 

PRINT  60aat(0ESCRPIII«I=l *281 
C PRINT  NUMBER  OE  FILTERS  USED 
PRINT  EQOltNUHFIL 

C PRINT  C8LCULATE0/RAW  FILTER  DAT* 

IF  INOCOFl  .EQ.  II  00  TO  10 
PRINT  BOlltNOCOFl 
GO  TO  20 
10  CONTINUE 

PRINT  6010  tNOCOFl 
20  IF  INUNFIL  .EQ.  II  GOTO  25 
IF  IN0C0F2  .EG.  II  GOTO  11 
PRINT  e013tN0C0F2 
GOTO  25 

11  PRINT  6012.N0C0F2 
C PRINT  FILTER  CALIBRATION  FROM  FILTAB 
25  IX=FILC*L*1 

PRINT  6040.  IFILTABCIX.IT I.IY=1. 71 
C PRINT  DAT*  TYPE 

IC*L:C*LD*T*1 
GO  TO  130.31. 33I.ICAL 
31  CONTINUE 

PRINT  6021.CALDAT 
PRINT  e023.C*LHDR 

6023  F0RMATIT20. 'CALIBRATION  : '.13*6.  A2  I 

GO  TO  40 
30  CONTINUE 

PRINT  6020.CAL0AT 
PRINT  6023.CALHDR 
GO  TO  40 
33  CONTINUE 

PRINT  G022.CALDAT 
40  CONTINUE 

C PRINT  GENERAL  WEIGHTING  TABLE 
IX:GENHEI*1 

PRINT  6050.(GFT*BIIX.IYI. IY=1.4I 
C PRINT  GENERAL  FUNCTION 

IFCGENFUN  .EG.  01  GO  TO  50 
PRINT  6031.GENFUN 
PRINT  6069.FUAVE.LMAVE 
PRINT  6070. OAT* 

GO  TO  60 
50  CONTINUE 

PRINT  6030.GENFUN 
60  CONTINUE 
C 

PRINT  6060.DFU.DFV 

6069  FORMATI T20. ’FOR  WAVELENGTHS  ».I4.'  TO  *.14  I 

6070  F0RMATIT20. 'GENERAL  FUNCTION  VALUE  IS  '.ES.?! 
RETURN 

C 

C «••••*••••••••••••••••*** 

c * * 

C * SUBROUTINE  TNITTL  * 


1 


B-23 


I 


• m 


I3ai* 

C 

• • 

1302* 

C 

1303* 

C 

1304* 

subroutine;  initil 

1305* 

C 

i3ae« 

C 

THIS  SUBROUTINE  INITIALIZES  NAIN  PROBRAN  VARIABLES 

1307* 

C 

1308* 

00  20  I = 1*2 

1308* 

00  10  J = lt340 

13ia» 

EINSTRIJ*!!  = 0 

1311* 

EIF0FXIJ.il  = 0 

1312* 

ELANBOIJ.II  = 0 

331  3» 

LLAHBDfJ.il  = 0 

1314* 

ERETLBIJ.II  = 0 

1315* 

GLAHBDIJ.il  : n 

1316* 

10 

CONTINUE 

131  7* 

20 

CONTINUE 

131  3* 

00  40  I = 1.4 

131  9» 

DO  30  J = 1.340 

1320* 

EFTIJ.II  = 0 

1321* 

30 

CONTINUE 

1322» 

40 

CONTINUE 

1323* 

DO  50  I = 1.340 

1324« 

ERETLBII.3I  = 0 

1325» 

DELTAdl  = 0 

1326* 

50 

CONTINUE 

1327* 

00  60  I = 1.9 

1328* 

FTLTER(I)  = 0 

1329* 

60 

CONTINUE 

1333* 

RETURN 

1331» 

C 

1332* 

c 

1333* 

c 

• # 

1334* 

c 

. SUBROUTINE  INTER”  * 

:335» 

c 

• • 

1336* 

c 

1337* 

c 

133S» 

SUBROUTINE  INTF RPf X . T , I. J .2 1 

1339* 

c 

134  9* 

c 

THIS  SUBROUTINE  INTEPOLES  THE  DATA  IN  ARRAY  T TO  CORRESPOND 

1341» 

c 

X AND  PUTS  THE  RESULT  IN  7 

1342* 

c 

1343* 

REAL  XI340.2I 

1344* 

REAL  Y(340.2I 

1345* 

A = YCJ.ll  - YIJ-l.ll 

1346* 

B : YIJ.2I  - VIJ-1.21 

1347* 

C = 3 / A 

1348* 

D r Xfl.ll  - YIJ-l.ll 

1349* 

” r 0 * C 

1350* 

; : YIJ-1.21  ♦ E 

1351* 

RETURN 

1352* 

c 

1353* 

c 

1354* 

c 

• # 

135  5* 

c 

* FUNCTION  SUH  * 

1356* 

c 

• % 

1357* 

c 

i 

) 


i 


i 


B.24 


1358* 

C 

1 

1359* 

REAL  FUNCTION  SUHtt.X.DI  ^ 

1360* 

C 

1361* 

C 

THIS  FUNCTION  SUNS  THE  PRODUCTS  OF  EACH  ELEMENT  OF  E *X  *0  ARRAYS 

1362* 

c 

1363* 

DIMENSION  EI340*2I* 

1364* 

1 XI340.2A* 

1365* 

2 013401 

136  6* 

K = 1 ) 

- 1367* 

SUM  =0  ^ 

' 1360* 

DO  50  I = l.MAXELM  1 

1369* 

IFIEII.l}  .LT.  01  GOTO  50  - 

1370* 

IFIE(I.l)  .LT.  Xll.lll  GOTO  50 

1371* 

10 

IFIEII.ll  .EG.  XlK.in  GOTO  30 

1372* 

IFIEII.ll  .LT.  XIK.III  GOTO  20 

1373* 

K = K * 1 

1374* 

GOTO  10 

1375* 

20 

CALL  INTERPIE. X.I.K.RI 

1376* 

GOTO  40 

1377* 

30 

R = X(K.2I 

1378* 

40 

SUM  = SUM  ♦ Ell. 21  * R * Dill 

1379* 

50 

CONTINUE 

1380* 

K = 1 

1381* 

DO  80  I = l.MAXELM 

1382* 

IFIEII.ll  .GE.  01  GOTO  80 

1383* 

DO  60  J = K.MAXELM 

1384* 

IFIEII.ll  .LT.  XIJ.lll  GOTO  70 

1385* 

60 

CONTINUE 

1386* 

70 

K = J 

1387* 

R = X(K.2I 

1388* 

IFIEII.ll  .NE.  XIK.III  CALL  INTE  IT>  IE  .X.  I . K HI  1 

1389* 

CONT  = Ell. 21  - E(I-1.21 

1390* 

IFKEII.ll  - Ell-l.lll  .EG.  IEII*1.1I  - Ell.  111!  CONT  = IEII*1.2I 

1391* 

1 ♦ E(I-1.2I 1 / 2 

1392* 

SUM  = SUM  ♦ CONT  * R • DII-ll 

1393* 

PK  = Ell. 21  - CONT 

1394* 

SUM  r SUM  ♦ PK  * R • Dill 

1395* 

80 

CONTINUE 

1396* 

RETURN 

1397* 

C 

1398* 

c 

1399* 

c 

1400* 

c 

* FUNCTION  SUMl  * 

1401* 

c 

• • 

1402* 

c 

1403* 

c 

mo*** 

1405* 

i4oe* 

1407* 

1408* 

1409* 

1413* 

1411* 

1412* 

1413* 

1414* 


REAL  FUNCTION  SUMl (E»DrL* H) 


-THIS  FUNCTION  SUMS  THE  PRODUCTS  CF  EACH  ELEMENT  OF  E tO  ARRAYS 


DIMENSION  Ef 340*21 t 
1 0(3401 

INTEGER  L*H 
SUMl  = 0 
DO  10  I = L*H 
IFIEII.ll  .LT.Ol  GOTO  10 
SUMl  : SUMl  ♦ Ed. 21  * D(II 


1415* 

10 

continue 

1416* 

DO  20  I = L.H 

I417* 

IFIEII.ll  .GT.  01  SOTO  20 

141S* 

CONT  : Ell. 2)  - EII-1.21 

1419» 

IFHEII.ll  - EII-1.1»I  .F8. 

lEII^l  d 1 

- FI  T,  H l» 

14Z0» 

1 ♦ E (I-1.2>l  / 2 

24ri» 

SUMl  r SUMl  ♦ CONT  • DII-ll 

i4rz» 

PK  : Ell, 21  - CONT 

14?3« 

SUMl  : SUMl  ♦ PK  » Dill 

14Z4* 

20 

CONTINUE 

i4r5» 

RETURN 

C 

l‘i?T 

C 

1428* 

C 

• 

• 

14rR» 

c 

• 

SU3R0UT  INE 

6LAMC0 

• 

I43a« 

c 

• 

• 

3 431  * 

c 

00  00 

143Z» 

c 

143  3* 

SUBROUTINE  GL AH  CO  1 R . X . Y . I 1 

1434* 

c 

1435* 

c 

SUBROUTINE  MULTIPLIES  EACH 

ELEMENT  IN 

X WITH 

ITS 

1436* 

c 

ELEMENT  IN  Y PUTTIN  IT  IN  ARRY  R 

1437* 

c 

1438* 

DIMENSION  R(340.2I. 

1438* 

1 XI340,2). 

144  3« 

1 Y(340.2) 

1441* 

K = 1 

1442* 

DO  50  I = l.MAXELM 

1443» 

10 

IFIXII.l)  .EQ.YlK.in  GOTO 

30 

1444* 

IFIXII.II  .LT.Y(K,1I)  CCTO 

20 

1445* 

K = K » 1 

144D» 

GOTO  10 

1447» 

20 

CALL  INTERPlX.Y.I.K.ZI 

1448» 

GOTO  40 

1443» 

30 

Z : YIK.2I 

1450» 

40 

R ( I. 11  : X( 1,1  ) 

14F1» 

RII, 2)  = X(T.2>  .2 

1452* 

50 

CONTINUE 

:4E  3» 

j:4 

14E4» 

RETURN  J 

1455* 

c 

• • •• 

1456* 

c 

• 

0 

1457. 

c 

• 

SU3R0UT  IN  E 

FILSIAI 

0 

1458* 

c 

• 

0 

1453. 

c 

» • •• • m 

00  00 

SUDRPUTINE  FILSUHI R.NtXtY  tYl.ZtD •$ I 

THIS  SU3R0UTINF  SUHS  THE  PRODUCTS  OF  EACH  1 EMENT  OF  XtY.YltZ.C  ARRAYS 


DIMENSION  Rn40tZ). 
1 XdtO.ZJt 

? Yt34ntZlt 

T YH340t2lt 

4 7l34Dt7tt 


1472» 

5 0(3401* 

1473* 

7 RLSI340*2I* 

1474» 

6 T(340*2t 

1475* 

DO  10  I = lf340 

1476* 

iFivti*ii  .re.  01  GOTO  20 

1477* 

Tdtll  = Tlltll 

1478* 

T(I<21  = VIIt2l  * V1(I.2I 

1479* 

10 

CONTINUE 

1480* 

20 

L = 8 

1481* 

SOTO  30 

1482* 

ENTRY  ETLSUBIRtNtX.T.ZtO.SI 

1483* 

L = 7 

1484* 

30 

M = N - 1 

1485* 

RINI  = 0 

148B* 

J = 1 

1487* 

K = 1 

1488* 

DO  110  I = ItHXXFLH 

1489* 

40 

IFIABSIXIItin  .ES.  TfK.lM  6070  6 0 

149Q* 

IFIABStXII.lll  .LT.  TIK.III  GOTO  50 

1491* 

K = K *1 

1492* 

GOTO  40 

1493* 

50 

CALL  INTERPlX.T.I.K.AI 

1494* 

GOTO  70 

1495* 

60 

A = T(K.2I 

1496* 

70 

IFCABSIXIItin  .EQ.  2(Jtin  GOTO  90 

1497* 

IFIABSCXII.lli  .LT.  2(J*11I  GOTO  80 

1498* 

J = J *1 

1499* 

GOTO  70 

isao* 

80 

CALL  INTERPfXtZtI. J*3I 

1501* 

GOTO  100 

1502* 

90 

B = ZIJ.2I 

1503* 

100 

RLSII.lt  = ABSfXII.lll 

1504* 

RLSd.ZI  = A * B 

1505* 

110 

CONTINUE 

1506* 

R(NI  = SUNCXfRLStOI 

1507* 

RETURN  L 

1508* 

C 

1509* 

C 

••  00000000  00  00  0 0 

1510* 

C 

0 • 

1511* 

C 

* FUNCTION  SEQUEN  * 

1512* 

C 

0 • 

1513* 

C 

1514* 

C 

1515* 

REAL  FUNCTION  SEQUEN 1 X .1. Yt J 1 

1516* 

C 

1517* 

c 

THIS  FUNCTION  RETURNS  THE  VALUE  IN  Y(I»2)  FOR  WHICH  YII.l)  HATCHES  Xd.ll 

1518* 

c 

1519* 

DIHENSION  X(340.2lt 

1520* 

1 T(340t2) 

1521* 

10 

IFCXII.ll  .EQ.  YfJ.in  GOTO  20 

1522* 

IF(X{I.1>  .LT.  Yej.lH  GOTO  30 

1523* 

J = J ♦ 1 

1524* 

GOTO  10 

1525* 

20 

SEQUEN  = Y(Jt2l 

1526* 

GOTO  40 

1527* 

30 

CALL  INTERPIXtYtl. J.SEQUENI 

1528* 

40 

RETURN 
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15?<9* 

c 

i5^a» 

c 

• ••••••••  m ••  •• 

1531» 

c 

• 

• 

:532» 

c 

• SUBROUTINE 

PRTCON  • 

1533* 

c 

• 

0 

153<*» 

c 

0 0 0 0 00  0 0 0 0 

1535* 

c 

1536* 

SUBROUTINE  PRTCON 

1537» 

c 

1638* 

c 

CONTROL  THE  PRINTING  OF  THE  LISTING 

1539* 

c 

1S40» 

LYNE  = LYNE  ♦ 2 

1541* 

IFCPEAK  .EC.  01  GOTO  20 

1542* 

PRINT  10 

1543* 

13 

format  I *0’  . 34 t •••!//  » 

1544* 

LYNE  3 LYNE  ♦ 4 

1545* 

PEAK  3 0 

1546* 

20 

IFILYNE  .LT.  LINMAXI  GOTO  30 

ir.4  7» 

PAGE  = PAGE  ♦ 1 

1548* 

CALL  HEAOIN 

:549« 

LYNE  = 12 

isr  3» 

33 

IFfEINSTRtT.H  .GT.  01  GOTO  60 

15E1* 

IFdLYNE  ♦ SI  .LT.  LINMAX  > GOTO  W 

1552* 

PAGE  = PACE  ♦ 1 

15t3» 

CALL  HEAOIN 

15E4* 

40 

PRINT  50 

1555* 

50 

FORMAT  C ‘O'/ZIOI •••1.5I*  •J.*PEAK».  51* 

*1.10!  *•*!  1 

15E6* 

LYNE  = LYNE  ♦ 4 

:E57» 

EINST3«I.ll  = A3S(EINSTR( r.l 1 t 

lEse* 

PEAK  : 1 

15E  9» 

60 

RETURN 

1563» 

c 

1561« 

c 

1562» 

c 

• 

0 

1563* 

c 

• SUBROUTINE 

HEAOIN  • 

15E4* 

c 

• 

• 

1563» 

c 

0 00  0 00  00  0 0 

1566* 

c 

15E7» 

SUEROUTINE  HEAOIN 

1SE8* 

c 

15&9» 

c 

PRINTS  THE  PAGE  HEADING 

157Q» 

c 

:57i» 

PRINT  10. DATE. PAGE 

:572» 

PRINT  lO.OrSCRP 

1573* 

PRINT  30. EVENT 

1574* 

PRINT  35. OMEGA 

1575* 

35 

FORMAT!*  SOURCE  SOLIC  ANGLE  IS  * .E  9. 21 

1576* 

PRINT  40 

I577* 

PRINT  50 

1578* 

PRINT  60.1 PRTl  AMfGENWEI»l .T) . 1=1 .3  » 

1579* 

PRINT  70 

1583* 

10 

FORMAT (* 1* .T2. *LMD  SPECTRAL  WEIGHTING 

PROG  » M*  .T90.*0A  TE  •.2AR. 

isei* 

ITllO.'PAGE  *.I3> 

15e2» 

20 

FORM  ATI *0* . T2.1 3A6. A2/T2. 13AG. A2 ) 

:583» 

30 

FORMAT!*  '.T2.3A6I 

lEe4» 

40 

FCRMAT!*0’.T31.  *AOJUSTEO  SPECTRAL 

SPECTRAL  SPECTRAL -RE  tin 

:585» 

lAL  GEWRAL  SPECTRAL  FILTER  T RA  NSNI SS  IONS*  1 
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158B* 

so 

F0RH8TI*  <iT3t*U«VE  CALIBRAT  INSTRUHENT  I NSTRUHENT 

1587* 

IRAOIANCE  IRRAOIANCE  UEIOHTING  FILTER 

1588* 

2 FILTER'! 

1588* 

60 

rORHATI*  'tTI^'LENSTH  FACTOR  READINOS  READINGS 

I530» 

1 SOURCE  (3-HNI  17-MHI  • t2  A6  tA?« Till)  »•  ONE 

1531* 

2NE  • TWO') 

1532* 

70 

FORHATI'  •tT2t' • .1 11 2X.3I •-» II 1 

1533* 

RETURN 

1S34* 

C 

1535* 

C 

1536« 

C 

1537* 

C 

• SUBROUTINE  BDREAO  • 

1538» 

C 

1533» 

C 

1600* 

C 

leoi* 

SUBROUTIIC  BDREADITABLEI 

1602» 

c 

1603* 

c 

SUBROUTINE  READS  BIOOECK  DATA  THAT  IS  ON  SEPARATE  CARD 

1604* 

c 

1605« 

DIHENSION  TABLE(34at2l 

IGQG* 

DO  20  I = 1.340 

1G07* 

READ  IS.10tERR=30MTABLr(I.JI.J=1.2  I 

IGOS* 

10 

F0RHATIF4.3.T3.E3.2I 

1633» 

20 

CONTINUE 

IGIO* 

GOTO  SO 

IGll* 

30 

READ  10.401  H 

1612* 

40 

F0RNATIA3I 

1G13» 

IFIH  .EG.  'END'I  GOTO  70 

1G14* 

SO 

PRINT  60 

1G15* 

GO 

FORMATI'O  BIOOECK  DID  NOT  END  CORRECTLY') 

IGIG* 

STOP 

1617» 

70 

RETURN 

iGia* 

c 

1619» 

c 

i62a» 

c 

• • 

1621* 

c 

• SUBROUTINE  HPLOT  • 

1G22* 

c 

• ♦ 

1G23* 

c 

1G24* 

c 

1625* 

SUBROUTINE  HPLOT 

1G2G* 

c 

1G27» 

c 

SUBROUTINE  CREATES  CALCOHP  PLOT  TAPE 

1G28* 

c 

1623* 

DIMENSION  XAXISfe85I.YAXISt685l .IBUFISOOO ) 

1633* 

INTEGER  EVENTI31 

1G31* 

CALL  PLOTSIIBUF. 5000.201 

1G32* 

CALL  PLOTIO. .-3G..-3) 

1G33* 

CALL  PLOTIO. .2. .-3) 

1634« 

CALL  SVHB0L10..0..0.21.1THAEHA.X24  7S .LMDSWP. 30. .171 

1635» 

CALL  PL0TI2..-0.8.-3I 

lG3e« 

NPLOT  = 0 

1637* 

4 

REAO fl5.10.ENDr993)  EVENT .ISIZE. IN CREH .FIG  .C H 

1G38* 

10 

FORMAT! 3A6.I4.I2.A2.AS) 

1G33* 

EVENT2llOI=EVENTfll 

1643* 

EVENT2I11I=EVENT(2I 

1G41* 

EVENT2(12)=EVENT(3I 

1G42* 

READI15.20)  XAXISIll 

FILTER 


SOURCE 

TWO 
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1643* 
1644* 
1645* 
1646* 
1647» 
1648* 
164  3* 
1653* 
1651* 
1652* 
1653* 
1654* 
1655* 
1656* 
1657* 
1658* 
1659* 
1663* 
1661* 
1662* 
1663* 
1664* 
1665* 
1666* 
166  7* 
1668* 
1669* 
1670* 
1671* 
1672* 
1673* 
1674* 
1675* 
1676* 
1677* 
1678* 
1679* 
1683* 
1681* 
1682* 
1683* 
1684* 
1685* 
1686* 
1687* 
1688* 
1689* 
1690* 
1691* 
1892* 
1693* 
1694* 
1695* 
1696* 
1697* 
1698* 
1699* 


20  FORM 
L = 2 
K:3 
lOCL 
00  3 
READ 
X AXI 
XAXI 
L = 

K = 

30  CONT 
XAXI 

rsrz 

00  5 
READ 
40  FORM 
TAXI 
50  CONT 
IF  « 
63  IF  ( 
CALL 
GOTO 
60  CALL 
CALL 
5 NPLO 
DRAW 
CALL 
CALL 
CALL 
CALL 
CALL 
LABE 
CALL 
CALL 
CALL 
DRAW 
CALL 
CALL 
00  5 
X r 
CALL 
CALL 
X = 
FPN 

55  CALL 
CALL 
SET 
XAXI 
XAXI 
Y-A 
ISIZ 
ISIZ 
IF  f 
CALL 
ISIZ 
YAXI 
CALL 


FIRST  PLOT 

NEW  STACK  OF  PLOTS 


AT IF5.31 


TA  : INCREM/2 
0 I=2fISIZE 
115*201  XAXISfLI 
SIL  »=XAXIS«LI-IDELTA 
SIKI=XAXISILI 
L*  2 
K ♦ 2 
INUE 

SILI  = XAXISIK-2)  * INCREH 

E r ISIZE*2 

0 I r l.ISIZEtZ 

115*401  YAXISlIt 

ATfE9.4l 

S«I*1»=YAXISCII 

INUE 

NPLOT  .EG.  0>  GOTO  5 9 FIRST  PLOT 

M00»NPL0T*4»  .EG.  01  GOTO  GO  9 NEW  STACK  OF  Pi 

PL0TI-n.79*7.1*-3) 

5 

PLOTt 10.15*-36. *-3» 

PL0TC3. .1.2*-3I 
T = NPLOT  ♦ 1 

BORDER*  HAVING  ESTASLISHED  LOWER  LEFT  CORNER 

PL0T«3. *8. *21 

PL0TI10.5*8.*2I 

PLOT< 10.5*0.*  21 

OLOTIO.  .0.*2» 

PLCT«a.79*1.3*-3»  a ORIGIN  OF  AXES 

L GRAPH 

SYMBOL (0.  *-0.9*0.  14*EVENT2*  0.  *7  2) 

SYMBOL  1 0.98 *-0.9 *0.14 *FIC*a.*  21 
SYMBOL < 6.1 6 *-0. 9*0.1 4 *CM*0. *5  1 
.NUMBER*LABEL  X-AXIS 
PLOTI 3. .0.*3» 

PL0H9.  *0.*2I 
5 1=1*11 
tl-ll  • 0.9 
PLOT! X*3. .3> 

PL0TIX*-Q.06*2» 

X - 0.37 
= fI*U  * 100. 

NUMBER  IX* -0.26  *3.  14  .FPN*  0.*  -1  » 

SYMBOL (3. 5 *-0.48 *0.14* 15HUAVELENGTH  I NM 1*0. *151 
PARAMETERS  FOR  LINE  GRAPH 
S(ISIZE*1J  = 200. 

SIIS1ZE*21  = 100.  * 10. /9. 

XIS  ROUTINES 

E JOCKEYED  TO  ELIMINATE  REDUMJANT  SCALING 
E = ISIZE/2 

MODINPLCT.ZI  .EG.  01  GOTO  61 
SCALEITAXIS*6.*ISIZE*2I  9 LINEAR 

E = ISIZC*2 

S(ISIZE*2I  = YAXISIISTZE*3 1 
AXISI3..0..30HSPECTRAL  IRRADIANCE  fW/CM  /NMI* 


1700» 

1701» 

1702* 

1703* 

1704* 

17Q5* 

1706* 

1707» 

1708» 

1709* 

1710* 

1711* 

1712* 

1713* 

1714* 

1715* 

1716* 

1717* 

1718* 

1719* 

1720* 

1721* 

1722* 

1723* 

1724* 

1725* 

1726* 

1727* 

1728* 

1729* 

1730* 

1731* 

1732* 

1733* 

1734* 

1735* 

1736* 

1737* 

1738* 

1739* 

1740* 

1741* 


A 3a*6.*9a.*VAXISIISIZC*ll*rAXTSt  tSIZC*2l • 

CALL  SVHBOLI-O.36f3.94*n. a7*lH2*  90 .* II 
CALL  LINriXAXIS*VAXIStI5TZL*ltD*OI 
OOTO  63 

61  CALL  SCAL0IVAXIS*6.tISIZC*21  A LOOARITIWIC 

ISIZE  = ISIZE*2 
C ORAU*LABEL  V-AXIS 

CALL  PL0Tta.*0.<3l 
CALL  PL0Tta.*6.*2) 

CALL  5VHBOL(-a.43*0.9*0.14>30HSPECTRAL  IRRAOZANCE  CW/CN  /NHI  * 
A 90.*30l 

CALL  SYHBOLI-0. 535*4. 43*0. a7tlH2 *9  0. *1  I 
C PLOT  BOTTOM  TICK  AND  NUMBER 

CALL  PL0TI3. *0.*3I 
CALL  PL0TC-O.12*O.»2» 

CALL  NUMBER! -0.4*0. *3. 14* 10. *0.* -1  I 
FPN  = ALOGIOIVAXISIISIZE*!)! 

CALL  NUMBER ! -0. 135 *0. 105*0. 07* rPN*  0. *-ll 
C PLOT  REMAINING  TICKS  AND  NUMBERS 

II  = 6.  * yAXISIISIZE*3»  a NUMBER  OF  CYCLES 

YY  = 0. 

DO  57  I = 1*11 
DO  56  J = 2*9 

X = YY  ♦ ALOGIOIJI/YAXIS! ISI2E*3  I 
CALL  PL0TID.*X*3» 

56  CALL  PL0T(-0.0e*X*2l 

YY  = I/YIAXISCISIZE*3I  9 INCH-ORDINATE  (T  END  OF  ITH  CYCLE 

CALL  PL0TI3. *YY*31 

CALL  PL0TI-0.12*YY*2I 

CALL  NUMBERI-a.4*YY*0.14.10.*0.*-l  > 

FPN  z FPN  ♦ 1. 

57  CALL  NUMBER!-0.135*YY*0.105*D.07 *FPN*n.*-l  ) 

YAXISIISIZE*2I  = YAXIS(ISIZE*3I 

CALL  LGLINf X AXIS* Y AXIS* IS  175*1*0.* 0* II 
GOTO  4 

C END  PLOT  ROUTINE 

999  CALL  PL0TI12..-36. *-3l 
CALL  PL0T(D.*O.5*-3l 
CALL  PLOTIO.*35.*2I 
CALL  PL0TI2.*-36.*999I 
RETURN 
END 


END  OF  compilation 
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K 

“x 

''x 

v/ 

Tx 

TAp,x 

l/C^ 

^JL 

.30-31 

• 00 

•on 

• 00 

•00 

• 00 

•ao 

•10*01 

701, 

.33-01 

• 00 

• 00 

• no 

•00 

• 00 

• 00 

•10*01 

202, 

.%0-01 

• 30 

• 00 

• 00 
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• 00 
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•10*01 

?33. 

•%%-01 

• 30 
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• 00 
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• on 

•00 

•10*01 

73%. 
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• 00 
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• 00 

• 00 

•10*01 

rns. 
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• 00 

• 00 

• 00 

• 00 

• 00 

• 00 
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?nT. 
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713. 
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?:&. 
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,30 
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.11*30 

,00 

• on 

• 09 

•00 

• 00 

•00 
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71%. 

.11*31 

,00 

• 00 

• 00 

•00 

• 00 

•no 

•10*01 

773. 
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.00 
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• 00 
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• 00 
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• 00 
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, 
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.00 
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• 00 
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.15*73 

.00 

.01 
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.10*01 
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• 00 
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• 00 
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• 30 
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.00 

•10*01 

7 73  . 

. l§  *73 

.30 
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.00 

•10*01 

7‘’1. 

. 19*93 

• 00 

• on 

• 00 

•00 

• 00 

• 00 

•10*01 

*•’7. 

.20*90 

• 00 

.00 

• 03 

• 00 

• 00 

• 00 

•10*01 

713. 

.?:*33 

.30 

• OP 

• 00 

•00 

• 00 

•on 

.10*01 

-3%. 
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.10 

• 30 

•00 
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• 00 
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773. 
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.00 
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• 00 

• 00 

•10*01 

7’&. 

.23*00 

.03 

.10 
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•no 

• 00 

• 00 

•10*01 

7’’7. 

.74*03 

.30 

.10 

• 10 

• 00 

• 00 

•00 

.10*01 

77*. 

.25*70 

.30 

.1'’ 

. 'n 

• 00 

• 00 

• 00 

.10*01 

77%. 

.76*30 

• 30 

.10 

.01 

• 00 

• 00 

• 00 

•10*01 

'*%3. 

.77*33 

.56*11 

.10 

• 01 

•no 

• 00 

•no 

.10*01 

‘*‘11. 

.29*73 

.56*31 

.1'' 

• m 
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• 00 

• 00 
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. 19*30 

.56*13 
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• 00 
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.10*01 
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.57*00 

•a" 

• n 

• 00 

• 00 

<00 

•10*01 
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.33*33 

.57*10 

•n 

.00 

• 00 

• 00 

• 00 

.10*01 

74S. 

.35*93 

•58*13 

.‘'7 

• on 

•no 

.01 

• 00 

.10*01 

746. 

. 36*-a 

.58*00 

• O'* 

.01 

.00 

• 30 

• 00 

•10*01 

7h7, 

.38*33 
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• on 

. m 

•no 

• 00 

• 00 

•10*01 

743. 
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.37 
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•00 

• 03 

•on 
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74%. 
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.45*99 
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.00 

• no 

•on 

• 00 

• 00 

•10*01 

?!•%. 

.63*93 
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Technical  Guide  - Hazard  Analysis  of  Broad-Band  Optical  Sources 

appendix  C 

OPERATING  INSTRUCTIONS 

The  following  data  will  not  change  from  run  to  run:  S-,  U-,  A-,  T-,  T-A-, 

C-A-,  V-,  V*-,  B-LAMBDA,  X-BAR  LAMBDA,  Y-BAR  LAMBDA,  Z-BAR  LAMBDA,  P-445, 

P-535,  and  P-575.  Other  data  vary  with  the  run;  the  data  decks  must  terminate 
with  the  proper  END  card  and  should  be  labeled  on  top.  The  run  request  must 
be  checked  to  be  certain  that  the  data  required  by  the  calculation  control  card 
are  supplied. 

There  must  be  two  source  description  cards,  even  1f  one  or  both  are  blank. 

The  cal  comp  plot  requires  wide  paper. 

The  same  calculation  and  calibration  data  may  be  used  for  several  sources. 
However,  changing  any  calculation  or  calibration  data  during  a run  can  be 
accomplished  by  using  'NEXT'  as  the  event  number  on  the  event  card.  Control 
then  branches  to  the  beginning  of  the  program,  and  expects  two  project 
description  cards,  etc,  as  1f  this  were  a different  run  request. 
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MESSAGES/HALTS 


Messages 

' (CALIBRATION  DATA  NOT  ENDED  CORRECTLY' 

'BIODECK  DOES  NOT  HAVE  AN  END  CARD' 

'FILTER  DATA  WERE  NOT  ENDED  CORRECTLY' 

'GENERAL  FUNCTION  DOES  NOT  END  CORRECTLY' 

'SPECTRAL  READINGS  WERE  NOT  ENDED  CORRECTLY' 

Causes-  No  END  card  for  data  section.  Or  more  than  340  cards  In  section. 
Effect  on  program-  Print  error  message  and  abort  run. 

Act1on-~Tnsert  END  card  at  end  of  appropriate  section;  If  END  card  Is 
present,  check  spelling.  Or  remove  extra  cards.  And  rerun. 

Message 

'WAVE  length  XXXX  of  SPECTRAL  IRRADIANCE  DOES  NOT  MATCH  WAVE  LENGTH 
OF  READING  FOR  CALIBRATION  FACTOR' 

Cause-  As  In  message. 

Effect  on  program-  Print  error  message  and  abort  run. 

Action-  Correct  errors  on  data  cards  and  rerun. 

Message 

■DIVISION  BY  ZERO  IN  CALIBRATION  FACTOR  SECTION  AT  WAVE  LENGTH  XXXX' 
Cause-  Datum  In  second  set  of  raw  calibration  data  Is  zero. 

Effect  on  program-  Print  error  message  and  abort  run. 

Action-  Correct  error  on  data  card  and  rerun. 

Message 

■calibration  FACTOR  ERROR-NUMBER  RAW  DATA  NOT  MATCHED' 

Cause-  Second  set  of  raw  calibration  data  has  different  number  of  Items 
than  first  set. 

Effect  on  program-  Print  error  message  and  abort. 

Action-  Insert  missing  data  and  rerun. 

Message 

■ATTEMRT  to  divide  A CALIBRATION  FACTOR  OF  ZERO  INTO  ADJUSTED 
INSTRUMENT  READING  AT  WAVELENGTH  XXXX' 

Cause-  First  set  of  raw  calibration  data  contains  a zero. 

Effect  on  program-  Print  error  message  and  abort. 

Act1oh-~Cliange  zero  to  valid  value  and  rerun. 

Message 

■ALL  PROCESSING  COMPLETED' 

Cause-  Normal  completion  of  program. 

Effect  on  program-  Normal  completion. 

Actlon-TTone  required. 
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INSTRUCTIONS  FOR  PILLINC.  OUT  PORfIS 

The  data  for  the  LMD  Spectral  Weiqhtinq  Program  must  be  entered  on  a number 
of  coding  forms  in  order  to  input  the  information  to  the  computer.  Explicit 
instructions  for  filling  out  each  form  are  provided  below: 

A.  Cover  Sheet.  These  data  are  used  to  identify  the  information  for  future 
reference  and  determine  the  kinds  of  input  data  to  the  Data  Processing 
personnel.  Item-by-item  instructions  are  given  bolov/. 

1.  Project  Number.  Pill  in  local  project  number  which  will  identify 
the  spectral  source  or  filter. 

2.  Source  Description.  Briefly  give  source  description  and/or  project 
name  (20  words  or  less).  Adequately  describe  source  for  future  reference. 

3.  Calibration  Deck  Used.  List  number  of  calibration  deck.  This  deck 
may  be  on  file  with  the  computer  center  or  the  calibration  factors  may  be 
submitted  with  the  rest  of  the  forms. 

4.  Calculation  Control  Card.  Enter  code  number  in  each  box  which 
describes  the  kinds  of  data  which  will  be  run.  Two  filters  with  known 
spectral  characteristics  may  be  submitted  with  the  source  data  if  desired. 
Normally,  instrument  readings  are  submitted  which  permit  the  computer  to 
calculate  the  spectral  irradiance,  E^,  by  dividing  the  reading  by  the 
calibration  factor i however,  spectral  irradiance  values  may  be  the  direct 
input  without  a calibration  factor  If  desired.  One  specific  biologic  function 
weighted  against  the  source  spectrum  may  be  listed  spectrally  If  desired. 

One  of  the  following  codes  should  be  Inserted  In  columns  4 and  5 of  this 
card  to  specify  this  function. 


00  — None 

01  — — Ultraviolet  Irradiance  According  to  ACGIH  Action  Spectra 

02  — Uj^  — Ultraviolet  Irradiance  According  to  CIE  Action  Spectra 

03  — Aj^  — Ultraviolet  Irradiance  According  to  ANSI  Action  Spectra 

04  — T — Transmission  of  the  Ocular  Media  x Ej^ 

05  — T*A  — Transmission  of  the  Ocular  Media  x Absorption  of  Retina  x 

06  — — Reciprocal  of  ANSI  MPE  Weighting  Factor 

07  — — Photopic  Spectral  Luminous  Efficiency  x E-y 

08  — V^'  — Scotopic  Spectral  Luminous  Efficiency  x E\ 

09  — B.  — Blue  Light  Hazard  Function  x Ej^ 

10  — — Spectral  Tristimulus  value  (red)  x E;^ 

11  — 7;^  — Spectral  Tristimulus  value  (green)  x E;^ 

12  — — Spectral  Tristimulus  value  (blue)  x Ex 

13  — P445  — Dartnall  Nomogram  Absorption  Coefficient  for  Blue  x Ex 

14  — P535  — Dartnall  Nomogram  Absorption  Coefficient  for  Green  x Ej^ 

15  — ^575  — Dartnall  Nomogram  Absorption  Coefficient  for  Red  x Ex 


Ex 
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5.  Number  of  Data  Decks  Submitted.  Enter  the  number  of  individual  | 

source  spectrums  which  are  being  submitted  at  this  time  under  this  cover  1 

sheet.  I 

•'  6.  Check  List.  Be  sure  all  listed  information  is  submitted  either  I 

on  or  with  the  cover  sheet.  I 

B.  Source  Data  Sheet.  This  form  is  the  first  form  used  for  recording  | 

spectral  data  by  wavelength  — Instrument  reading  from  a particular  instru-  | 

ment.  All  three  items  on  this  form  must  be  filled  in.  J 

1.  Source  Name.  The  source  name  is  used  to  distinguish  between  several  ! 

source  spectrums  entered  under  the  same  cover  sheet.  Hiis  name  should  ; 

be  kept  under  8 characters  if  possible.  j 

j 

2.  Source  Solid  Angle.  The  Solid  Angle  of  the  source  must  be  entered  ! 

in  exponential  notation  (example,  l.OOE-05).  If  the  solid  angle  is  not 

known,  an  approximation  must  be  used.  j 

i 

3.  Data  Values.  Values  are  entered  in  integer  form  for  wavelength 
and  exponential  form  for  instrument  reading.  The  wavelength  values  must 
be  placed  in  the  rightmost  columns  ( , 235, not  ^35  , ).  Spectral  peaks 

may  be  identified  to  the  right  of  the  instrument  readings  by  inserting  > 

the  word  "PEAK"  in  the  indicated  columns . These  values  are  then  treated  . 

separately  by  the  computer  program,  ; 

C.  Calibration-Filter  Sheet.  This  form  may  be  used  for  calibration  or  for 
filter  information.  Information  should  be  placed  in  the  first  two  columns 
by  wavelength  and  value.  The  type  of  coded  information  should  be  identified 
at  the  top  of  the  form. 

D.  Data  Continuation  Sheet.  This  form  may  be  used  for  continuation  of 
either  source  data  or  calibration  values. 


E.  End  Card.  At  the  end  of  either  the  source  data  deck  or  at  the  end  of 
the  calibration  deck,  the  statement  "END  of  Spectral  Data"  or  "END  of 
Calibration  Deck"  must  be  inserted. 
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EXPANDED  PROGRAM  USE 

Some  additional  features  are  built  into  the  program.  However  their  use 
is  not  reconsnended  for  general  use.  Key  personnel  may  be  required  to 
assemble  the  run  stream  and  thereby  increase  the  processing  time.  Some  of 
these  features  are  listed  below. 

A.  Distance  Factors.  These  two  factors  may  be  used  to  adjust  calibration 
values  for  variations  in  the  distance  from  a standard  lamp  that  a piece 

of  equipment  was  placed  for  calibration.  For  instance,  if  the  required 
measurement  dist^ulce  '.’as  50  cm  and  the  calibration  was  made  at  100  cm, 
a distance  factor  of  4.0  could  be  inserted  into  the  program.  However, 
the  READ  statement  was  deleted  from  the  progr2un  which  read  in  these 
values.  These  two  factors  are  designated  DFU  and  DFV  and  are  applied 
to  the  UV  portion  of  the  spectrum  and  the  visible  and  near-IR  portion 
respectively. 

B.  Uncomputed  Calibration  Deck.  The  calibration  values  may  be  submitted 
in  two  decks  rather  thzm  one.  The  program  will  then  divide  the  instru- 
ment readings  (first  deck)  by  the  spectrum  of  the  standard  lamp  (second 
deck).  If  data  are  submitted  in  this  form,  a zero  is  placed  in  column  3 
of  the  Calculation  Control  Card. 

C.  Uncomputed  Filter  Data.  Filter  data  may  l)e  submitted  in  two  columns 
rather  than  one  following  the  wavelength  values.  The  second  column  is 
then  divided  by  the  first.  Enter  a "2"  in  column  2 of  the  Calculation 
Control  Card  to  run  the  program  this  way. 

D.  General  Function.  Part  of  the  source  data  may  be  corrected  if  found 
faulty  after  the  cards  have  Ijeen  punched.  The  deck  may  be  corrected 
throughout  certain  wavelength  regions  by  the  use  of  this  function.  The 
starting  wavelength,  the  ending  wavelength,  and  the  correction  factor 
are  coded  on  the  special  form  provided.  As  tmny  correction  cards  as 
needed  may  be  used.  An  "END"  card  must  be  placed  at  the  end  of  this 
function. 

E.  Expanded  Run  Stream.  The  order  of  the  various  inputs  is  listed  below. 

1.  Source  Description  (two  cards) 

2.  Calculation  Control  Card 

3.  Distance  Factors  (program  change  required) 

4.  Calibration  Deck 

5.  Standard  Lamp  Data  (if  required) 

6.  BIO  Deck 

7.  Filter  Transmissions 

8.  General  Function 

9.  Source  Name  (Event) 

10.  Source  Solid  Angle 

11.  Source  Data  Readings 
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LASER  MICROWAVE  DIVISION  SPECTRAL  WEIGHTING  PROGRAM 


COVER  SHEET 

A.  PROJECT  NUMBER;  42; • 

B.  SOURCE  DESCRIPTION:  


C. 


D. 


(free  format, 

158  spaces) 


.CALIBRATION  DECK  USED;  DECK  NUMBER 

□ ON  PILE  □ INCLUDED  WITH  DATA  DECK 

CALCULATION  CONTROL  CARD:  1234  5 


BOX  NUMBER:  I I I I I I 

1.  Number  of  Filters  to  be  Processed  — Enter  Number (maximum  is  two) 

2 . Leave  Blan)c 

3.  Form  of  Calibration  Dec)c; 

a.  Regular  — Enter  "1" 

b.  None  Required  — Enter  "2" 

4-5,  Specific  Biologic  Function  to  be  Listed  Spectrally  — Ejiter  code 


a 


a. 

None 

— 00 

e. 

T — 04 

i. 

^x' 

' — 08 

m. 

z — 

X 

12 

b. 

sx  — 

01 

f. 

T-A  — 05 

j. 

®X 

— 09 

n. 

P445 

— 13 

c. 

“x  -- 

02 

g. 

1/C^  — 06 

k. 

— 10 

o. 

^535 

— 14 

d. 

*x  -- 

03 

h. 

Vx  — 07 

1. 

^x 

— 11 

P- 

P575 

— 15 

•fUMKII  nr  DATA  DECKS  SUBMITTED:  . 

•e  » • r*r  f 1 Bmiros  'iMcrlptlan  Q Instruwant  Raadinga 

mr4  I ] An^l* 
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